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INTRODUCTION 


The  determination  of  physical  properties  of  liners,  explosives, 
and  propellants  is  important  for  the  characterization  of  cook-off 
behavior,  aging  processes,  and  the  mathematic  modelling  of  the  thermal 
responses  of  warheads  and  rocket  motors. 

The  development  of  test  hardware,  data  acquisition,  and  numerical 
methods  for  data  reduction  and  analysis  is  of  considerable  importance. 
Thermal  analysis  techniques  such  as  Differential  Scanning  Calorimetry 
(DSC),  Thermogravimetric  Analysis  (TGA),  and  Differential  Thermal 
Analysis  (DTA)  have  traditionally  been  used  to  determine  the  thermal 
characteristics  and  decomposition  kinetics  of  reactive  substances. 

As  compared  with  other  analytical  techniques,  thermal  analysis 
has  until  recently  received  little  attention  in  chemical  research. 
This  is  primarily  due  to  instrument  insensitivity  and  the  laborious 
task  of  manuaHy  reducing  data.  More  sensitive  techniques,  usually 
spectroscopic  in  nature,  were  opted  for  in  the  determination  of  chemi¬ 
cal  kinetics. 

There  has  recently  been  a  resurgence  of  interest  in  thermal 
analysis  due  to  the  increased  safety  awareness  of  both  the  industrial 
and  military  sectors.  This  is  due  to  the  development  of  new  materials 
for  use  in  energy  development,  and  for  extraterrestrial  and  military 
applications. 

Considerable  activity  has  been  generated  in  the  area  of  the 
application  of  known  numerical  techniques  to  thermal  analysis.  This 
has  been  made  possible  with  the  development  of  new  thermal  analysis 
techniques,  the  availability  of  microprocessor-controlled  analysis 
instrumentation,  and  the  development  of  small  computer  technology  for 
the  acquisition  and  reduction  of  large  quantities  of  data. 

As  a  part  of  an  ongoing  process  of  upgrading  the  Naval  Weapons 
Center's  analytical  capabilities,  the  Explosives  Formulation  Branch 
has  written  computer  software  to  acquire  and  analyze  data  from  labora¬ 
tory  instruments.  These  programs  include  routines  to  input  digital 
data  from  interfaced  digitizers  and  various  numerical  analysis  and 
graphics  subroutines  to  reduce,  display,  and  analyze  the  data. 

Chemical  kinetics  subroutines  have  been  developed  for  nonisother- 
mal  DSC,  isothermal  DSC,  and  for  TGA.  This  has  been  accomplished 
through  the  digitization  of  the  analog  outputs  of  the  laboratory 
instrumentation  which  is  achieved  by  interfacing  with  a  Fluke  7600A 
scanning  digital  multimeter  or  Fluke  8520A  digital  multiswter.  Trans¬ 
fer  to  the  Hewlett-Packard  9845T  computer  is  achieved  via  the  IEEE-488 
interface  bus  either  by  direct  computer  command  or  by  use  of  service 
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request  interrupts.  The  rate  of  data  acquisition  is  programed  via 
coanand  string  to  the  microprocessor  control  units  of  the  digital 
multimeters. 


THEORETICAL  BACKGROUND 


Mathematical  interpretation  of  thermodynamic  data  is  accomplished 
through  the  application  of  the  Gibbs-Helmholtz  equation,  usually 
stated  as: 

.  AH 


AC/T 

(1/T) 


where  G  is  the  Gibbs  free  energy,  H  is  the  enthalpy,  and  T  the  abso¬ 
lute  temperature.  Kinetic  data  may  be  determined  from  the  relation¬ 
ship  AG°  =  -RT  In  K(T),  K(T)  being  the  equilibrium  constant  at  temper¬ 
ature  T.  For  a  reaction  such  as  A  t  B,  the  equilibrium  constant  may  be 
written  K(T)  =  k2(T)/k2*(T)  where  k2(T)  is  the  rate  of  the  forward 
reaction  and  k2*(T)  is  the  rate  of  the  reverse  reaction. 

Arrhenius  first  gave  the  expression  for  these  rate  constants  as: 
k2(T)  -  Z  •  axp(-E*/*T) 

where  Z  is  an  experimentally  determined  rate  constant  with  the  dimen¬ 
sions  mole/sec  and  E*  is  the  activation  energy.  The  form  of  this 
equation  was  taken  as  a  logical  extension  of  the  Gibbs-Helmholtz 
equation  and  is  used  presently  in  experimental  determinations.  Cur¬ 
rent  interpretation  of  E*  is  that  of  an  energy  barrier  or  depth  of  a 
potential  well,  and  the  factor  exp(-E*/RT)  is  the  probability  per 
unit  time  that  a  molecule  will  penetrate  the  barrier. 

From  the  Gibbs-Helmholtz  equation,  the  temperature  dependence  of 
the  equilibrium  constant  is  given  as: 


Unfortunately,  when  dealing  with  highly  energetic  materials,  the 
chemical  reactions  are  irreversible.  The  temperatures  and  pressures 
are  not  uniform  and,  in  fact,  have  severe  gradients  throughout  the 
material,  presenting  a  problem  in  nonequilibrium  thermodynamics.  When 
dealing  with  small  sables,  temperature  and  pressure  gradients  are 
usually  assuattd  to  be  negligible  and  a  rate  equation  may  be  written 
for  the  reaction.  For  the  reaction  B  Products: 

Ute  -  $  -  k(T)  •  f(X) 
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where  X  is  the  fraction  of  B  that  has  reacted  and  f(X)  is  some  func¬ 
tion  which  depends  upon  the  mechanism  by  which  the  reaction  takes 
place. 

For  simple  liquid  and  gaseous  reactions,  f(X)  takes  the  form 
(1-X)n,  where  n  is  referred  to  as  the  reaction  order,  n  is  usually  2 
or  less.  A  set  of  coupled  reactions  will  yield  fractional  n's  while 
reactions  which  occur  in  the  solid  state  via  diffusion  or  nucleation 
mechanisms  take  on  complex  algebraic  forms  with  n  usually  less  than  1. 

When  looking  at  energetic  materials,  n=l  works  quite  well  in  the 
higher  temperature  ranges  where  B-»  Products  is  an  accurate  representa¬ 
tion  of  the  reactions  taking  place.  However,  at  lower  temperatures, 
nucleation  reactions  are  present.  They  are  dependent  upon  the  surface 
free  energy  of  the  material  which  is  affected  by  such  factors  as 
physical  state,  impurities,  and  method  of  preparation.  Thus,  a  study 
of  reaction  rates  at  higher  temperatures  may  give  a  somewhat  limited 
view  of  the  manner  in  which  a  material  might  behave. 

For  nonhomogeneous  materials  such  as  propellants,  reaction  rates 
are  diffusion  controlled  and  highly  pressure-dependent.  For  materials 
of  current  interest,  f(X)  will  stand  for  (1-X)  for  simple  reactions 
and  X(l-X)  for  the  autocatalytic  reactions  which  occur  in  many  explo¬ 
sives. 

For  materials  which  melt  with  decomposition,  complex  kinetics 
occur  in  the  initial  stages  of  the  decomposition  reaction.  Behavior 
beyond  the  exotherm  maximum  is  usually  more  typical  of  the  first  order 
reaction,  ®( liquid)-*  Products.  The  first  stages  of  the  decomposition 

process  frequently  draw  the  most  interest  in  the  thermal  stability 
studies.  Investigation  of  this  behavior  will  give  a  more  complete 
description  of  how  a  particular  substance  reacts. 


ANALYTICAL  TECHNIQUES 


GENERAL 

Ordinate  data  points  are  obtained  by  digitizing  millivolt  output 
from  the  instrumentation  with  the  sensitivity  of  the  instruments 
selectable  by  the  operator.  The  analog  output  of  the  DSC  represents 
energy  absorption  (endotherm)  or  release  (exotherm)  in  units  of  milli- 
calories  per  second.  The  TGA  output  represents  milligrams. 

Abscissa  data  points  are  obtained  as  time  from  the  Real  Time 
Clock  interfaced  to  the  HP  9845T  computer.  Temperature  is  obtained  as 
the  product  of  time  and  heating  rate.  Extremely  accurate  temperature 
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tracking  has  been  achieved  in  the  case  of  the  DSC,  while  less  favor¬ 
able  results  have  been  obtained  from  the  TGA.  This  is  due  to  the 
heating  rates  deviating  from  programmed  values.  The  problea  has  been 
eliminated  by  manually  entering  data  into  a  separate  computer  program 
after  the  TGA  run  is  complete. 

Data  points  (Xj,  Yj)  are  curve-fitted  using  a  precise  continuous 
cubic  polynominal,  or  cubic  spline,  algorithm.  Ordinate  values  are 
represented  as: 


where  the  coefficient  matrix  M  is  determined  by  demanding  continuity 
of  the  polynominal  and  its  first  derivative  at  each  data  point. 

Use  of  this  algorithm  allows  facile  differentiation  and  integra¬ 
tion  of  the  data  curves  and  is  relatively  insensitive  to  noise  in  the 
instrumental  output. 

NONISOTHERMAL  DSC 

Graphics  output  of  a  DSC  run  of  RDX  Type  II  Class  5  is  shown  in 
Figure  1.  Figure  2  is  the  data  after  a  linear  baseline  correction  has 
been  made. 

For  a  nonisothermal  DSC  run,  the  ordinant  displacement  dH/dt  is 
given  as: 

D(t)  -  AH  •  rata  -  AH  4—  •  AH  •  k(t)(1  -  X)"  (1) 


where  AH  is  the  total  enthalpy  change  for  the  reaction,  represented  as 
the  area  A,  under  the  displacement  -  temperature  curve.  Thus  dX/dt 
=  1/A  •  dH/dT.  If  this  equation  is  integrated,  the  result  is  X 
=  AHx(T)/A  where  AHx  is  the  heat  released  due  to  the  fraction  X  react¬ 
ing.  This  is  represented  by  the  area  a(T),  under  the  peak  up  until 
the  temperature  in  question.  Thus  X(T)  =  a(T)/A  may  be  determined 
graphically  or  by  a  numerical  integration  technique  as  is  done  by  the 
computer  program  which  has  been  developed  (Figure  3). 

Borchardt  and  Daniels1  have  shown  that  the  rate  constant  may  be 
expressed  as: 

_  Ml)  - 

'Hans  J.  Brochardt  and  Farrington  Daniels.  "The  Application  of 
Differential  Thermal  Analysis  to  the  Study  of  Reaction  Kinetics," 
Journal  of  the  American  Chemical  Society,  Vol.  79,  No.  4,  July  1956, 
p.  41 . 
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Comparing  the  rate  constants  at  two  temperatures  gives: 

k(T2)  (dH/dt) (T2 )  |A-a(T,  )]  " 

k?V  "  (dH/dt)  (T, )  *  [a-.(T2)J 

Since  dH/dT  is  the  ordinate  deflection  D(T),  and  X  =  a(T)/A,  this 
relation  may  be  rewritten: 

k(Ta)  D(Ta)  p-x(Tt)l" 

k?V  “  D?V  ‘  [l  -  x(T2)J 


Assuming  only  one  chemical  reaction  takes  place, 


k(T)  may  be  written: 


k(T)  -  Z  •  exp(-E*/RT) 

and  thus: 


ft*  /J_  ,x1  »<V 

•»»(-«* rtiJT  "  “'Ll  ' T,  ’  Ti  'J  "  id,) 


Thus  the  activation  energy  may  be  calculated  from  ordinate  dis¬ 
placement  and  partial  areas.  In  fact,  for  closely  spaced  data  points, 
the  last  term  approaches  zero.2 

The  frequency  factor  Z,  and  thus  the  rate  constant  may  be  calcu¬ 
lated  using  partial  areas  and  known  values.  With  $  as  the  heating 
rate : 

-  Z  •  exp  [-I*/ET(r)J  [1  -  *&-] 


In  the  computer  implementation  of  this  method,  E*  is  calculated 
from  a  linear  least  squares  (linear  regression)  algorithm  which  fits 
the  curve  of  In  D  versus  1/T  (Figure  4). 


2 Jen  Chiu.  "Dynamic  Differential  Calorimetric  Technique  for 
Measuring  Heat  of  Polymerization,"  in  Proceedings  of  the  Second 
Symposium  on  Analytical  Calorimetry,  ed .  by  Roger  S.  Porter,  New  York 
Plenum,  1970,  p.  171-183. 
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Figure  5  shows  a  plot  of  the  reaction  rate  versus  fraction  decom¬ 
posed.  Ideally,  the  plot  should  provide  some  idea  of  the  algebraic 
form  of  the  rate  equation,  with  a  straight  line  indicating  a  first 
order  reaction.  For  an  autocatalytic  reaction,  a  straight  line  would 
result  from  plotting  rate  versus  X(l-X).  In  the  case  of  solids,  there 
are  several  processes  occurring  simultaneously  and  a  clear  and  precise 
interpretation  of  this  plot  would  be  difficult  to  achieve.  In  a  homo¬ 
geneous  liquid,  such  as  is  present  to  the  right  of  the  exotherm  maxi¬ 
mum  in  the  case  of  RDX,  a  portion  may  be  represented  by 

RDX  (Liquid)"*  Pro<*ucts: 

lata  -  k(T)  (1  -  X)  (2) 


and  appears  as  a  virtually  linear  portion  of  the  curve.  This  portion 
may  be  selected  for  determination  of  kinetics  for  the  first  order 
reaction  of  Equation  2. 

Kissinger3  has  devised  a  method  based  upon  variable  heating  rates 
to  determine  the  activation  energy.  This  method  is  independent  of 
reaction  order  and  should  give  an  accurate  value  for  E*  assuming  that 
only  one  chemical  reaction  is  taking  place  at  a  time.  The  basic 
equation  is  given  as: 


d(la  ♦/T«)/d(l/T)  -  -E*/R 

where  Tm  is  the  temperature  at  which  the  maximum  exotherm  or  endothera 
occurs.  A  plot  of  In  $/Tm  versus  l/Tm  from  several  runs  will  have 
slope  -E*/R. 


ISOTHERMAL  DSC 


The  software  developed  for  nonisothermal  DSC  has  been  modified  to 
allow  for  the  treatment  of  the  isothermal  case. 


Because  the  temperature  is  constant  in  an  isothermal  DSC,  the 
equation: 


dX 

dt 


11*  «(t2)  -  *(tt) 

ta-e,  A  •  (t,  -  t,) 

•*o 


may  be  solved  directly  for  k(T). 


k(T) 


f(x) 


(3) 


3 Homer  E.  Kissinger.  "Reaction  Kinetics  in  Differential  Thermal 
Analysis,"  Analytical  Chemistry,  Vol.  29,  June  1957,  p.  1702. 


8 


*■**«»  ^wwpx 


NWC  TP  6337 

For  closely  spaced  data  points,  dX/dt  is  reasonably  approximated 
as  Aa/At*A.  A  plot  of  Aa/A/At  versus  a/A  will  indicate  *he  nature  of 
f(X),  and  k  is  calculated  directly  by  linear  regression  analysis  for 
the  linear  portion  of  this  curve  when  first  order  kinetics  are 
observed.  One  could  plot  rate  versus  different  algebraic  forms  of 
f(X)  to  see  if  a  straight  line  results.  Again,  f(X)  would  be  expected 
to  take  the  form  ( 1 -X)  or  X(l-X)  for  explosives. t  Figures  6  to  9  are 
graphic  output  examples  of  this  program  and  Tables  1  and  2  provide 
examples  of  numeric  outputs. 

The  activation  energy  is  calculated  by  plotting  the  values  of 
k(T)  obtained  at  several  temperatures  versus  1/T.  The  slope  of  this 
plot  is  E*/R. 

Previous  isothermal  kinetic  determinations  have  used  the  portion 
of  data  to  the  right  of  the  exothermic  peak.  By  integrating  Equation 
3,  a  plot  of  fraction  unreacted  versus  time  should  have  the  form: 

1  -  X  ■  exp  [-k(T)  •  tj'  (A) 

for  a  first  order  reaction  (Figure  10). 4 

From  Equation  1,  the  ordinate  displacement  D(t)  is  given  as 
AH  *  k(T)  •  (1-x).  For  a  first  order  reaction: 

D(t)  •  AH  4“f  •  AH  •  k(T)  •  (1  -  *) 


Substituting  Equation  A  for  (1-x)  gives: 

D(t)  -  AH  •  k(T)  •  expt-k(T)  •  t] 

or  ln(D(t)l  -  In [AH  •  k(T)]  -  k(T)  •  t 


tjamey  and  Rodgers.  Los  Alamos  Scientific  Laboratory,  Private 
Communication. 

“Ray  N.  Rogers.  "Differential  Scanning  Calorimetric  Determination 
of  Kinetics  Constants  of  Systems  That  Melt  With  Decomposition,"  Thermo- 
chemica  ACTA,  Vol.  3,  1972,  p.  437-A47. 
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Thus  a  plot  of  la  D  versus  tine  should  have  slope  k(T)  for  a 
first  order  reaction.5  Figure  11  shows  the  plot  of  In  D  versus  time 
for  RDX  Class  5.  The  linear  portion  of  this  curve  was  used  for  deter¬ 
mination  of  first  order  kinetics  for  RDXq^-»  Products. 

From  the  fora  of  Equation  1,  it  is  seen  that  In  D  =  n  ln(l-a) 
+  constants,  so  a  plot  of  In  D  versus  ln(l-a)  should  yield  the  reac¬ 
tion  order  n.  Figure  12  represents  this  plot.  Calculation  of  the 
slope  in  the  region  from  which  kinetic  data  were  calculated  gives  a 
value  of  1.05,  which  agrees  with  the  assumption  of  a  first  order  reac¬ 
tion. 


Both  the  partial  area  technique  and  the  old  technique  of  plotting 
In  D  versus  t  give  essentially  the  same  answer  in  the  region  chosen. 
The  partial  area  technique  is  more  dynamic  in  that  the  fraction  of 
material  reacted  can  be  calculated  at  any  point.  Thus,  if  some  reac¬ 
tion  mechanism  is  hypothesized,  it  may  be  verified  with  actual  experi¬ 
mental  numbers. 

THERMOGRAVIMETRIC  ANALYSIS 

While  the  DSC  data  are  useful  for  describing  the  thermodynamics 
of  a  system,  TGA  is  necessary  for  correlating  weight  loss  with  tem¬ 
perature.  For  reactions  of  the  type  solid  -*•  gas,  TGA  may  be  used  for 
kinetic  determinations.  The  fraction  X,  given  as  X  =  M-m/M,  M 
=  original  mass,  and  the  reaction  rate  dx/dt  nay  be  taken  directly 
from  the  TGA  data. 

Numerical  analysis  for  TGA  data  is  the  same  as  that  for  noniso- 
thermal  DSC.  Because  the  temperature  rise  rate  in  the  TGA  instrument 
is  nonlinear,  temperature  must  be  calculated  after  the  run  is  com¬ 
pleted.  Temperature/mass  data  are  then  manually  entered  into  the 
computer.  Temperature  is  accurately  determined  by  entering  millivolt 
values  from  strip  chart  output  and  converting  to  temperature  values 
using  a  standard  National  Bureau  of  Standards  interpolating  polyno¬ 
mial.  Data  for  RDX  Class  5  is  shown  in  Table  3  and  is  plotted  in 
Figures  13  through  15. 

This  program  was  developed  primarily  to  investigate  the  decom¬ 
position  of  rocket  motor  liners  where  time,  temperature,  and  mass 
values  are  important.  If  only  mass  and  temperature  are  needed,  the 
software  is  easily  modified  to  allow  entry  of  approximately  ten  times 
the  number  of  data  points.  The  acquisition  of  computer  compatible 
thermogravimetric  equipment  will  also  allow  acquisition  of  more  data 
points  when  all  three  reference  points  are  required. 


5 Ray  N.  Rogers.  ’'Simplified  Determination  of  Rate  Constants  by 
Scanning  Calorimetry",  Analytical  Chemistry,  Vol.  44,  December  1971, 
p.  1336. 
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CONCLUSIONS 


The  acquisition  and  analysis  of  large  amounts  of  data  with  a 
digital  computer  has  opened  a  new  field  in  the  area  of  thermal  anal¬ 
ysis.  Plotting  and  numerical  analysis  software  allow  the  data  to  be 
reduced  in  a  fraction  of  the  time  previously  needed  and  permit  use  of 
a  larger  number  of  data  points. 

In  addition,  fast  integrating  algorithms  allow  use  of  partial 
area  techniques  which  will  prove  to  be  very  influential  in  future 
analytical  procedures. 

The  successful  application  of  these  computer  techniques  in  the 
determination  of  isothermal  kinetics  should  prove  extremely  useful  in 
determining  the  cook-off  hazards  from  aging  studies  and  understanding 
cook-off  behavior. 
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FIGURE  1.  DSC  Plot  of  RDX  Class  5.  (Scanning  Rata  is  10  dag/ain.) 
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FIGURE  2.  Plot  of  RDX  Class  5  with  a  Llnaar  Baseline  Correction 
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FIGURE  7.  Isothermal  Run  of  RDX  Class  5  op  Perkin  Elmer  DSC  II.  (The  shaded 
area  represents  the  area  calculated  by  the  computer.  Endpoints  for  the  deter¬ 
mination  as  well  as  points  used  for  baseline  determination  are  operator  selec¬ 
table.) 


FIGURE  8.  Plot  of  the  Rate  of  Change  of  the  Area  Fraction  vs  the  Area  Fraction. 
(The  linear  portion  of  the  curve  Is  chosen  for  the  determination  of  kinetic 
parameters . ) 
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I.  Graphics  of  Theraogravlaetrlc  Analysis  of  SOX  Class  5.  (Scan, 
dag/aln.) 
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TABLE  t .  Computer  Printout  of  Dot*  Used  to  Determine 
Kinetic  Psrsmeters.  (The  section  of  the  data  that 
is  used  for  the  calculations  is  operator  selec¬ 
table.) 


RD;;  Hus  5  10  dtg-nin  2  «c*l'»4c  2.32  »g 

OctoOtr  2,  1981 
KinttKi  ton 


T«>c  *>  1  T  Art*  Fr  tc  t . 


LN<ordt  rt*t  t  >  R«1  >t  i  R»M  Pt> 


479 

2. O9E-03 

0.004 

480 

2.08E-O3 

0.003 

481 

2.086-03 

8.006 

482 

2. 086-03 

0.007 

482 

2.076-03 

0.803 

48  i 

2.076-03 

0.010 

484 

2.076-03 

0.011 

483 

2-066-03 

0.013 

485 

2. 066-03 

0.813 

438 

2.066-03 

0.017 

487 

2.056-03 

0.020 

487 

2.036-03 

0.022 

433 

2.03E-03 

O.02S 

439 

2.056-03 

0.028 

430 

2.046-03 

0.032 

430 

2.046-03 

0.033 

491 

2.046-03 

0.039 

492 

2.036-03 

0.044 

493 

2.036-03 

0.049 

493 

2.036-93 

0.034 

494 

2.026-03 

0.060 

495 

2.026-03 

0.066 

498 

2.026-03 

0.073 

496 

2.01E-03 

0.080 

497 

2.016-03 

0.088 

498 

2.016-03 

0.097 

498 

2.016-03 

0.  106 

49? 

2.006-03 

0.116 

500 

2.006-03 

0.  127 

501 

2.00E-03 

0.139 

581 

1.996-03 

0.  132 

502 

1.996-03 

0.163 

5  03 

1.996-03 

0. 179 

504 

1.996-03 

0. 194 

504 

1.986-03 

0.210 

505 

1.986-03 

0.227 

506 

1.986-03 

0.243 

507 

1.976-03 

0.264 

307 

1.976-03 

0.284 

508 

1.976-03 

0.305 

509 

1.976-03 

0.327 

50? 

1.9C6-03 

8.  349 

510 

1.986-83 

0.373 

311 

1.986-03 

0.390 

312 

1.936-03 

0.423 

312 

1.936-03 

0.449 

513 

1.936-03 

0.475 

314 

1.936-03 

0.301 

315 

1.946-03 

0.527 

315 

1.946-03 

0.334 

SIC 

1.946-03 

0.380 

317 

1.946-03 

0.606 

518 

1.936-83 

0.631 

318 

1. J36-03 

8.636 

31? 

t.936-83 

0.660 

.001 

1.306-03 

5 

.001 

1.306-03 

6 

.126 

1.366-03 

7 

.283 

1.57E-03 

8 

.415 

1 . 806-03 

9 

.323 

2.026-03 

10 

.646 

2.24E-03 

11 

.734 

2.496-03 

12 

.830 

2.71E-03 

13 

.923 

2.976-03 

14 

1.029 

3.276-03 

IS 

1.117 

3.606-03 

16 

1.203 

3.91E-03 

17 

1.279 

4.246-03 

18 

1.372 

4.39E-03 

19 

1.461 

3.03E-03 

20 

1.557 

5.30E-03 

21 

1.647 

6.046-03 

22 

1.743 

6.616-03 

23 

1.840 

7.286-03 

24 

1.922 

7.956-03 

25 

2.001 

8.606-03 

26 

2.080 

9.316-03 

27 

2. 161 

1.016-02 

28 

2.242 

1.09E-O2 

29 

2.318 

1.186-02 

30 

2.393 

1.276-02 

31 

2.481 

1 . 386-02 

32 

2.549 

1.496-02 

33 

2.624 

1.606-02 

34 

2.693 

1 . 726-82 

33 

2.734 

1.836-02 

36 

2.814 

1.956-02 

37 

2.874 

2.076-02 

38 

2.932 

2.196-02 

39 

2.986 

2.326-02 

40 

3.039 

2.446-02 

41 

3.092 

2.386-02 

42 

3.142 

2.71E-02 

43 

3. 180 

2.836-02 

44 

3.223 

2,976-02 

45 

3.269 

3.086-02 

46 

3.312 

3.236-02 

47 

3.347 

3. 36E-02 

48 

3.372 

3.466-02 

49 

3.368 

3.536-02 

50 

3.307 

3.566-02 

51 

3,397 

3.566-02 

52 

3.405 

3.606-02 

33 

3.399 

3.616-02 

54 

3.309 

3.566-02 

55 

3.374 

3.536-02 

56 

3.356 

3.476-02 

57 

3.332 

3.406-02 

56 

3.296 

3.306-02 

39 
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TABLE  1.  (Coned.) 


320 

1.92E-03 

0.  784 

3.261 

3. 13E-02 

60 

520 

1.92E-83 

0.726 

3.220 

3.076-02 

61 

521 

1 . 92E-03 

0.748 

3. 173 

2. 93E-02 

62 

522 

1.92E-03 

0.768 

3.130 

2.306-02 

63 

523 

1.91E-93 

0.  738 

3.033 

2.63E-02 

64 

323 

1.  916-03 

0.806 

3.034 

2.3 6E-02 

63 

524 

1.91E-83 

0.824 

7.992 

2.446-02 

66 

525  • 

1.91E-03 

0.841 

2.932 

2.  34E-02 

67 

526 

1.90E-03 

0.  833 

2.921 

2.266-02 

68 

526 

1 . 90E-03 

0.374 

2.899 

2. 20E-O2 

69 

527 

1.90E-03 

0.890 

2.383 

2.  16E-02 

70 

328 

1.89E-03 

0.906 

2.  375 

2. 146-02 

71 

329 

1.89E-03 

0.921 

2.867 

2.  12E-02 

72 

529 

1 . 89E-03 

0.936 

2.849 

2. 09E-O2 

73 

530 

1.896-03 

0.931 

2.803 

2.02E-02 

74 

531 

1.88E-03 

0.  963 

2.673 

1.87E-02 

75 

331 

1.88E-03 

0.976 

2.  376 

1.33E-02 

76 

332 

1.88E-03 

0.984 

1.920 

1 . 03E-02 

77 

333 

1.88E-03 

0.989 

1.443 

6.47E-03 

78 

334 

1.87E-03 

0.992 

.990 

4.07E-03 

79 

534 

1.87E-03 

0.994 

.313 

2.38E-03 

80 

333 

1.87E-03 

0.995 

.032 

1.6OE-03 

81 

336 

1.37E-03 

0.996 

.001 

1. 18E-03 

32 

337 

1.366-03 

0.996 

.001 

1.22E-03 

83 

337 

1.866-83 

0.997 

.001 

1.2OE-03 

84 

338 

1.86E-03 

0.998 

.001 

1.20E-03 

83 

339 

1.86E-03 

0.999 

.001 

1 . 206-03 

36 

339 

1 . 33E-03 

1.000 

.001 

1.206-03 

87 

28 
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TABLE  2.  Computer  Printout  of  Data  Uaad  for  Graphlca 
and  Kinatlc  Paraaatar  Calculations. 


TIME  »  • 

AREA 

FRACT  A 

AREA  INCREMENT 

RATE  OF  INCREASE 

Ft  8 

95 

800 

8. 133 

213.982 

23.067 

3 

1«4 

1007 

a.  193 

206.693 

23.963 

6 

113 

1207 

0.233 

200.093 

23. 106 

7 

121 

1398 

0.272 

191.348 

22.124 

8 

138 

1381 

0.38? 

182.243 

21.05? 

9 

1  39 

1734 

8.341 

17  5.670 

20.039 

10 

14? 

1919 

0.373 

134. 933 

19.061 

11 

133 

2873 

8.484 

156.694 

18.094 

12 

135 

2223 

0.432 

148.303 

17. 191 

13 

173 

2337 

0.  430 

142.030 

16.403 

14 

132 

2502 

0.483 

134.818 

13.372 

13 

191 

23.9 

0.311 

127.723 

14.734 

16 

193 

2  750 

0.335 

120.343 

13.939 

17 

208 

2334 

0,35? 

114.329 

13.202 

10 

217 

2973 

0.370 

100.207 

12.499 

19 

223 

3073 

0.398 

102.329 

1 1.819 

20 

234 

3172 

0.613 

93.783 

11.183 

21 

242 

3234 

0.334 

91.899 

10.614 

22 

23 1 

3351 

0. 331 

87.433 

10.108 

23 

230 

3434 

0.338 

83.27? 

9.617 

24 

238 

3313 

8.383 

78.900 

9.  123 

23 

27? 

3389 

0.390 

75.347 

8.702 

26 

233 

3330 

0.711 

71.571 

8.269 

27 

294 

3728 

0.723 

38.069 

7.061 

20 

303 

3793 

0.737 

33.031 

7.312 

29 

312 

3853 

0.749 

32. 191 

7.  183 

30 

320 

3915 

0.731 

59.341 

6.87? 

31 

329 

3972 

0.772 

37.043 

6.387 

32 

338 

4027 

0.783 

54.783 

6.329 

33 

343 

4080 

0.793 

52. 794 

3.099 

34 

335 

4130 

0.803 

30.711 

3.836 

33 

334 

4179 

0.912 

43.098 

3.646 

36 

372 

422? 

0.822 

47.532 

3.491 

3? 

381 

4273 

0.831 

45.991 

3.313 

38 

390 

4317 

0.839 

44.413 

3.  129 

39 

398 

4331 

8.648 

43.472 

5.021 

40 

407 

4*0? 

0.833 

42,417 

4.397 

41 

413 

4445 

0.834 

41.302 

4.795 

*2 

424 

448.5 

0.872 

40.514 

4.678 

43 

433 

45;?3 

0.879 

39.628 

4.377 

44 

442 

45t.4 

0.897 

30.837 

4.486 

43 

430 

*402 

0.894 

38.838 

4.393 

46 

439 

4339 

0.982 

37.534 

4.338 

47 

438 

4373 

0.909 

36.961 

4.269 

40 

478 

4712 

0.916 

36.206 

4.  191 

49 

485 

4748 

0.923 

33.731 

4.  129 

30 

494 

4784 

0.930 

33.279 

4.074 

31 

302 

4818 

0.937 

34.702 

4.808 

32 

311 

4832 

0.943 

34.209 

3.931 

S3 

320 

4883 

0.930 

33.339 

3.876 

34 

328 

4919 

0.  933 

32.772 

3.786 

33 

33? 

4931 

0.932 

31.932 

3.690 

36 

343 

4982 

0.939 

30. 791 

3.556 

37 

334 

3811 

9.974 

29. 172 

3.368 

30 

333 

3038 

0.979 

26.901 

3.107 

39 

371 

3831 

0.984 

£3.672 

2.734 

60 

380 

3881 

0.986 

19.826 

2.290 

61 

389 

3897 

0.991 

16.039 

1.833 

62 

597 

3118 

9.993 

12.610 

1.437 

63 

303 

3128 

0.993 

9.003 

1.132 

64 

313 

3127 

8.997 

7.634 

0.002 

63 

323 

3131 

0.990 

3.936 

9.606 

66 

332 

3133 

0.999 

4.307 

9.530 

6? 

341 

3141 

0.999 

3.329 

8.303 

69 

349 

3143 

1.990 

2.  106 

0.243 

69 

338 

3144 

1.098 

1.339 

9.193 

70 

337 

3143 

l.AfC 

0.346 

9.040 

71 

29 
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TABLE  2.  (Contd.) 


111 

Ti««« i«c  > 

V 

Loo(Y) 

4 

<l-«> 

La  100»> 1-* 

s 

93 

24.4 

3.193 

0. 136 

0.044 

4.436 

6 

104 

23.3 

3.138 

0. 196 

0.804 

4.387 

7 

113 

22.6 

3.  120 

0.233 

0.763 

4.336 

8 

121 

21.6 

3.071 

0.272 

0.723 

4.288 

9 

130 

20.6 

3.023 

0.307 

0.693 

4.238 

10 

139 

19.6 

2.973 

0.341 

0.639 

4.  138 

11 

147 

18.6 

2.922 

0.373 

0.627 

4.138 

12 

136 

17.6 

2.868 

0.  404 

0.596 

4.088 

13 

163 

16.0 

2.822 

0.432 

0.368 

4.039 

14 

173 

16.0 

2.771 

0.460 

0.340 

3.989 

1? 

182 

13.2 

2.719 

0.406 

0.314 

3.939 

16 

191 

14.3 

2.662 

0.311 

0.409 

3.090 

17 

199 

13.6 

2.607 

0.333 

0.463 

3.040 

18 

200 

12.0 

2.333 

0.337 

0.443 

3.792 

10 

217 

12.2 

2.498 

0.370 

0.422 

3.743 

20 

223 

11.5 

2.441 

0.390 

8.482 

3.693 

21 

234 

10.9 

2.380 

0.616 

0.384 

3.647 

22 

242 

10.3 

2.336 

0.634 

0.366 

3.399 

23 

231 

9.9 

2.289 

0.631 

0.349 

3.331 

24 

260 

9.4 

2.236 

0.660 

0.332 

3.304 

29 

269 

8.9 

2. 108 

0.603 

0.317 

3.437 

26 

277 

0.3 

2.  138 

0.698 

0.302 

3.409 

27 

286 

0.1 

2.086 

0.711 

0.289 

3.362 

29 

234 

7.7 

2.039 

0.723 

0.273 

3.313 

29 

303 

7.3 

1.994 

0.737 

0.263 

3.268 

30 

312 

7.0 

1.930 

0.749 

0.231 

3.221 

31 

320 

6.7 

1.906 

0.761 

0.239 

3.174 

32 

329 

6.4 

1.864 

0.772 

0.228 

3.126 

33 

338 

6.2 

1.627 

0.783 

0.217 

3.079 

34 

346 

6.0 

1.788 

0.793 

0.207 

3.030 

35 

353 

3.7 

1.747 

0.803 

0.197 

2.981 

36 

364 

3.6 

1.717 

0.812 

0.  188 

2.932 

37 

372 

3.4 

1.639 

0.822 

0.178 

2.802 

38 

381 

3.2 

1.649 

0.831 

0.169 

2.830 

39 

390 

3.1 

1.626 

0.839 

0.161 

2.778 

40 

398 

5.0 

1.608 

0.848 

0.  132 

2.724 

41 

407 

4.8 

1.379 

0.836 

0. 144 

2.660 

42 

416 

4.7 

1.535 

0.864 

0. 136 

2.611 

43 

424 

4.6 

1.531 

0.872 

0. 128 

2.331 

44 

433 

4.5 

1.512 

8.879 

0. 121 

2.489 

43 

442 

4.4 

1.489 

0.887 

0.113 

2.423 

46 

430 

4.4 

1.474 

0.894 

0. 106 

2.337 

47 

439 

4.3 

1.460 

0.902 

0.098 

2.205 

48 

468 

4.2 

1.442 

0.909 

0.091 

2.209 

49 

476 

4.2 

1.423 

0.916 

0.084 

2.129 

30 

405 

4.  1 

1.412 

0.923 

0.077 

2.042 

31 

494 

4.0 

1.396 

0.930 

0.070 

1.949 

32 

502 

4.8 

1.382 

8.937 

0.063 

1.840 

33 

311 

3.9 

1.363 

0.943 

0.037 

1.737 

34 

520 

3.8 

1.343 

0.930 

e.ese 

1.616 

33 

328 

3.7 

1.320 

0.936 

8.044 

1.400 

36 

337 

3.6 

1.290 

0.962 

0.030 

1.320 

37 

546 

3.5 

1.244 

0.968 

0.032 

1.133 

38 

334 

3.3 

1.181 

0.974 

0.026 

0.930 

39 

563 

2.9 

1.878 

0.979 

0.021 

0.734 

60 

371 

2.3 

0.922 

0.984 

0.016 

0.483 

61 

590 

2.1 

0.726 

0.908 

8.012 

0.214 

62 

389 

1.6 

0.497 

8.991 

0.009 

-0.076 

63 

397 

1.3 

0.240 

0.993 

0.007 

-0.303 

64 

686 

1.0 

-0.003 

0.993 

0.003 

-0.711 

65 

615 

.0 

-0.253 

0.997 

0.003 

-i.oro 

66 

623 

.6 

-0.508 

0.990 

0.002 

-1.400 

67 

632 

,3 

-0.760 

0.999 

0.001 

-1.970 

60 

641 

.3 

-1.213 

0.999 

0.001 

-2.600 

69 

649 

.2 

-1.367 

1.080 

0.008 

-3.419 

70 

659 

.1 

-2.416 

1.000 

0.000 

-3.003 
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TABLE  3.  Computer  Printout  of  Time,  Temperature 
Lose  Data  for  TGA  Program. 


*Q‘<  Cl  is*  5  TGA  3  dt$'i»in  . 

•  1*1  «ui  is  21.39  H**tmg  ~it*  is  3  d*g^*m 

•#n*»t»vM»y  is  10  ng/'iiv 


'ttpirituf*  <C' 

T  i  *+<*> i  n) 

Hats(iftC) 

nut  Let* 

36 

0.00 

21.6 

0.00 

37 

2.50 

21.6 

0.00 

38 

5.00 

21.6 

0.00 

41 

7.38 

21.6 

0.00 

45 

10.00 

21.6 

0.00 

49 

12.58 

21.6 

0.00 

33 

13.00 

21.6 

0.00 

61 

17.30 

21. * 

0.00 

€7 

20.00 

21.6 

0.00 

74 

22.30 

21.6 

0.00 

82 

25.00 

21.6 

0.00 

90 

27.30 

21.6 

0.00 

99 

30.00 

21.6 

0.00 

108 

32.50 

21.6 

0.00 

117 

33.00 

21.6 

0.00 

126 

37.50 

21.6 

0.00 

134 

40.00 

21.6 

0.00 

1*2 

42.30 

21.6 

0.00 

149 

45.00 

21.6 

0.00 

137 

47.30 

21.6 

0.00 

154 

30.00 

21.6 

0.00 

171 

52.30 

21.6 

e.00 

179 

35.00 

21.6 

0.00 

137 

57.30 

21.6 

0.00 

194 

60.00 

21.6 

0.80 

201 

62.30 

21.2 

1.85 

208 

bS .  00 

20.3 

3.09 

214 

67.30 

17.8 

17.60 

224 

70.00 

10.4 

51.88 

231 

72.30 

2.9 

86.61 

239 

73.00 

.4 

98.  19 

246 

77.30 

.  t 

99.38 

236 

80.00 

0.0 

188.80 

264 

82.50 

0.0 

100.00 

272 

83.90 

0.0 

100.00 

281 

87.30 

0.0 

100.00 

291 

90.00 

0.0 

100.00 

Mass 


Data  Pt  t 


1 

2 

3 

4 

3 

6 

7 

e 

9 

le 

it 

12 

13 

14 
13 
16 
17 

15 

19 

20 
21 
22 

23 

24 
23 
26 

27 

28 

29 

30 

31 

32 

33 

34 
33 

36 

37 
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PROGRAM  FOR  AGUISITION  AND  ANALYSIS  OF  THERMOGRAV I  METRIC  BATA 


18  OPTION  BASE  1 

28  COM  T<  300;,Mas*c3O0).M*««l<  300) , Ti m*<  300  • ,  N,  M8.  ft*  C  88  1 ,  B*  1 50 1 
25  COM  A' 300 > , Bc380> , An«u< 300) ,Bn*w<300> , D*r i v(300> 

38  COM  As'0:9)  i  Co«ffici«nts  For  int*rpol at i no  polynomial 
40  MPT  READ  Aa1  for  t ypt  S  thtrmoc oupl *. 

45  LINPUT  “ENTER  RUN  IDENTIFIER", BS 

50  P 1 ot  *  1 

60  IF  P1o»>1  THEN  LINPUT  "BO  YOU  WANT  TO  PLOT  ANOTHER  RUN?",Cont* 

62  IF  iFlo»/P  AND  <Cont SC  1 , 13< >"Y“ >  THEN  GOSUB  Cont  l  nut 
65  INPUT  "ENTER  NUMBER  OF  DATA  POINTS", N 

70  REDIM  T<N),Ma*sCN>,Nassl<N>,Tia*\N) 

88  INPUT  "ENTER  MILLIVOLT  READING  FOR  INDIUM  MELT " , Mp 

98  INPUT  “ENTER  THERMOCOUPLE  CAIN  (1  OR  5>“,Cain 

100  INPUT  “ENTER  INITIAL  MASS‘,M0 

105  INPUT  "ENTER  HEATING  RATE" , Hrat* 

110  INPUT  "ENTER  MILLIVOLT  VALUE  FOR  INITIAL  MASS" , Mo 

115  INPUT  "ENTER  MASS  SENSITIVITY  Cmgmv):  1  or  10  ",S*ns 

120  PRINTER  IS  16 

130  FOR  !«t  TO  N 

140  PRINT  LINC6) 

150  PRINT  USING  220; 1 

160  INPUT  X,Y 

170  X»<X-Mp)'G*i n*l .  882  I  Calculat*  S  thsrmocoupl*  output, 

180  T< I >«FNSeal <X' 1000)  !  Calculat*  T*mp*ratur«. 

190  Nass<I>«M0*S*nt*<Y-Mv)  1  Calculate  Itui  For  c Sens  >*grmv  sensitivity. 
195  Ti *»< I )■< I-I )*2. 5  1  Tim*  in  minut**:  5  minrinch  :hart  sp**d. 

200  Massl < T )«CM0-Mass< I ' ) “M0*I00  <P«rc*nt  Mass  Lot*. 

202  PRINT  PAGE 

205  PRINT  USING  225; TCI ) , Hasstl ) , T i »•(  I > 

210  NEXT  I 

220  IMAGE  "ENTER  MILLIVOLT  READINGS  FOR  POINT  »",BB,”  C T«»p. , Hass ) “ 

225  IMAGE  "T*»o*ratur«  is  *,DDD,"  C;  Mass  is  ",DD.DD,"  mg.;  Tim*  is  ‘ 
. D , "  min." 

230  GOSUB  Printer 

235  LINPUT  “DO  YOU  WANT  THE  DATA  STORED?" , St or«S 

236  IF  Stor«*Cl,  M*"Y"  THEN  CALL  Stor*<TC* > , Mass< « ■ , M0, N,  BS) 

240  IF  Plot>!  THEN  GOSUB  S 

245  REM  Plot  X  Mass  loss  vs  Ttmpsratur* 

250  INPUT  "ENTER  PLOTTER  CODE:  1«CRT;  2* 9872A", Plotter 
260  IF  Plott*r*l  THEN  PLOTTER  IS  "GRAPHICS" 

270  IF  Plot* *r»2  THEN  PLOTTER  IS  7,5, "9872A" 

280  GRAPHICS 

290  LOCATE  15,90,20,70 

295  INPUT  "ENTER  HIGHEST  PERCENTAGE  DESIRED  ON  PLOT”,Q 

300  SCALE  0,600,0,0 

310  LINPUT  “ENTER  PLOT  LABEL", AS 

320  Xnam*S»" TEMPERATURE  <C>" 

3 30  Yna**S»"P*re*nt  Msss  Loss* 

340  AXES  100, QyiO, 0, 0 

350  FRAME 

368  MOVE  Til), Massif  1 > 

370  FOR  1-2  TO  N  1PL0T  FIRST  ARRAY. 

300  DRAM  TCI), Mass  1 < I ) 

390  riEXT  I 

400  PENUP 

410  LORG  6  •  Lab*!  X  Axis 

420  FOR  t>0  TO  600  STEP  100 
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430  MOVE  1,-2 

440  LABEL  I 

450  NEXT  I 

466  LORG  3  1  Late)  V 

470  FOR  1=0  10  TO  u  STEP  0-10 

430  MOVE  0,  I 

430  LABEL  I 

500  NEXT  I 

510  lORG  1  '  Lib*  I  Plot. 

510  MOVE  0,1.03*0 

530  LABEL  A* 

540  LORG  6 

541  MOVE  300, -.03*0 

542  LABEL  XriimiS 

543  MOVE  25,0-2 

544  LDIR  PI-2 

545  LABEL  Y nan it 

54?  PENUP 

543  LDIR  O 

550  PAUSE 

553  IF  P1o»t*=*l  THEN  DUMP  GRAPHICS 

555  LINPUT  "WOULD  YOU  LIKE  TO  REPEAT  THIS  PLOT-” . Rp* t 
56u  IF  Rpt*Cl,lJ»"Y*  THEN  GOTO  250 
565  plot -Plot* I 

563  GOTO  60 

570  i:  1  Plot  oth*r  arrav-s 
580  IF  P 1 ot *2  THEN  LINE  TYPE  4 

530  IF  Plot =3  THEN  LINE  TYPE  6 

600  IF  P I ot  *4  THEN  LINE  TiPE  8 

610  MOVE  Tvl>,Ma*sl< 1) 

620  FOR  I «2  TO  N 

630  DRAW  T' 

640  NEXT  1 

650  PENUP 

660  p I ot  =P 1 ot ♦ l 

670  GOTO  60 

703  E-it:  IF  P’ott*r=2  THEN  OUTPUT  ?,5;"SP0-IH" 

710  IF  Flott«r*2  THEN  STOP 
720  DUMP  GRAPHICS 

730  EXIT  GRAPHICS 

740  GCLEAR 

7' 9  END 

760  Prinuri  PRINTER  IS  O 
765  PRINT  PAGE 
770  PRINT  LIN'2' 

’80  PRINT  B* 

735  PRINT  USING  790; MO, Hr  at #, Sins 

790  IMAGE  "Initial  miss  is  ",DDS.SD,"  1*3.  Ht  at  i  ng  r  it  *  is  ".DD,"  d«g  mo",  , 
"Sins  1 1 1  hi  t  y  is  ",  DD,'*  »g-»v" 

300  PRINT  L I N v  2 i 

310  PRINT  *  Ttroptr  it  urt  ''Cl  T  i  mi  ( w  i  n  )  Miss  (  mg  )  \  Mass  Loss 

Dai  j  Pt  »"  "  ' 

320  PRINT 

330  FOP  I«1  TO  N 

335  IP  Miis< I '0  THEN  Miss<I>«0 

336  IF  Missl1  l  /  >  1 00  THEN  Mass  1  •  I  ■ » 100 

340  PRINT  USING  360; TC I >, Ti»«- I >,Mass < I I.Miss 1 < I >, I 
850  NEXT  I 

360  IMAGE  DDDD!6:<,DDD.DDSX,DD.D10X,DDD.DBI2X,DDD 

370  PRINT  PAGE 

380  RETURN 

1 000  PATA  .927763167, 1 69526. 51 5, -  31 56336  3. 94 , 39907 306 33, - 1 . 63S65E I  2 . 1.83027E14 
.  -1 .  77241  E  16. 6.  1750!  El  7, -1.56  1 05-El  9,  l .  695  35E20 
1065  DEF  FliSc *1  < 2 1 
!0’5  COM  Ai0:9> 
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10S5  M*#rpol  At.  ing  pel  ynmi  il  for  converting  milliv-olt*  to  d«gr*«i 
1095  •  or  typ*  S  tpl it i nu»-10Lrhodt urn  plji  mu»  t  her  no:  oup  1  *  . 

1105  :  :  ii  =A<:0  >*Z*<  Ac  1  »S*<At2  >*2*tA'  3  t*II»(.  At  4  )*2*tRv.5>  *2*<rt<6  7  8 ■> 

Z*Rv9  ■  '»,>>> 

1115  PETURN  So il 
1125  -  --EHD 

1155  :  .6  St  or* ' X< * • , Yt * ) , Y0, N, Id* > 

1 165  t-OPT  A'  N  >.E  <  M) 

1175  n-T  A*:< 

1135  "AT  B  =  Y 

1135  .INPUT  "ENTEP  FILE  NAME  " ,  D* 

1210  .INPUT  "ENTER  MASS  STORAGE  MEDIUM", E* 

1215  :»£AT£  B*« "4£*, N  16»2 

1225  ASSIGN  #4  TO  D*8."  1  "&E* 

1235  'PINT  »4;ld*,Y0,N 
1245  MAT  PRINT  *4;  A  ,  E 
1255  ASSIGN  »4  TO  * 

1235  SJBEXIT 
1275  SUBEND 
1300  Id  Continue 

1305  ;pt i on  ease  i 

1310  ISM  T  -  *>,«»**>*>,«***  l<*>,Ti  »«■  4»,N,  MO,  ft*.  S* 

1320  ICM  A 1  ♦  ■ ,  E ( *  ■' ,  An*  m  c  *  1 ,  Bn*w  *1  *  >  ,  Otr  1  o  (  *> 

1325  7IM  P »g»* C 30 1 

1  370  IALL  PI  ot  t  *r  <B*.  Ttm*t:*>,Mi*»l  <  *  > ,  N,  2> 

1340  PRINT  PAGE 

1350  .INPUT  “BO  YOU  WANT  TO  CONTINUE  WITH  DATA  ANALYSIS?" , V* 

1360  IF  V * C 1 , 1 3 <> " Y "  THEN  SUBEXIT 

1370  I-vPUT  "ENTEP  BEGINNING  AND  ENDING  POINT  FOR  DATA  ANALYS  IS "  ,  P 1 ,  P2 
1380  I  t*,«P2-Pl  +  l 

1  330  *  ED  I M  A<  Dt  m>  ,  B>  Di  ,  D«r  i  ut  Di  »-l  )  ,  An«w<Di  m- 1  >  ,  Bn««<Di  »-l  > 

1400  FOR  K*  1  TO  Dim 
1410  a.K>»TCR1+K-1> 

1420  E ■ K )»M**j  npi »K-I ' 

1430  v.£i<T  K 

1440  PjP  1*1  TO  Dim-1  1  Cilculit*  dtriuitiuis  it  midpoint*. 

14  50  Anew V I  > * <; A '  I ♦  1  )+A(  1  > > r2 
1460  -,EXT  I 

1470  Epj*lE-6 
1480  =PINT  PACE 
1 480  PRINT  L  I  N  C  3  > 

1500  PRINT  B* 

1510  PRINT  "Sub**t  Output" 

1520  PRINT  L  I N 1 2  > 

15  30  IALL  Spl  1  n*<  Dt  m,  Di  m-1 ,  AOO  ,  B*  *)  ,  An*w<*  )  ,  Bnei-i*:  *  )  ,  D*r  i  y  <  *  )  ,  0 1 ,  Eps  > 

1540  "'AT  D»r  i  u»RBS>  D«r  i  u) 

1550  R  id*  *»”  T«mp«r  it  ur«  (C  1  Miis  <  md )  -dM-  dT  1  nCdM^dT  **  * 

1560  =RINT  Pig«* 

1570  PRINT  L I N< 2 1 
1580  FOR  1*1  TO  Dim-1 

1590  fpinT  USING  1 670;  I ,  An«w<  1 )  ,  Bn«u<  I  >  ,  D«r  1  v<  I  > ,  LOG<  D*r  i  v<  I  :■  > 

1600  IF  I  MOD  5000  THEN  GOTO  1660 
1610  PRINT  PAGE 
1628  =RINT  Bt 
1630  PRINT  L IN t  2 ) 

1640  PRINT  Pig** 

1650  PRINT  L IN ( 2 ) 

1660  NEXT  I 

1670  IMAGE  DDDXXX, DDDDXXXXXX, MDD.  DDDXXXXXX,  MD.  DDEXXXXXX,MZ.  DDD 
1630  IALL  Plott*r<B*,A<*>,B<*),Dim,0> 

1690  CALL  Plotter'  8* , An*v( # I ,  D«r  1  ut  *!> ,  Di  m-1 ,  11 
1700  'PINT  PAGE 
1710  PRINT  LIN<2.' 

1720  .INPUT  "DO  YOU  WANT  TO  CALCULATE  k«"ETIC  PAPAMETEPS  PROM  THE  DATA?",S* 
1730  IF  S*CI,1J*"Y"  THEN  CALL  n i n«t i c»( E* , finivf  • ' , Bn«u' *1 , D*r lu<*> , Di a, M0> 
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1740  GOTO  1356 
1750  SUBEND 

ISOO  SUB  Plot t «r< A*, X« *> , Y(»> ,N.Fl4g > 

1801  INPUT  "ENTER  PLOTTER  CODE:  1*CRT  2*98?2A" , Pl ot t «r 
1302  IF  Plo*.t«r*2  THEN  PLOTTER  18  ?,5,"9872A" 

1803  IF  FI »g*2  THEN  G**“Ti  »,*■  mi  nut  «i  >  " 

1305  IF  F 1 4g< 2  THEN  G**"T*»p«r»tur«  <C>“ 

1806  IF  Fl*g*0  THEN  H«*>>ni  11  igrmt* 

1307  IF  F1»g*l  THEN  H**"M i I • i gr  ioj  d« g“ 

1808  IF  F!»g*2  THEN  H*»"P*i~i:«nt  M4*i  Lo**“ 

1310  PRINTER  IS  0 

1811  PRINT  PAGE 

1812  PRINT 

1813  PRINT  A* 

1814  IF  Fl*g*l  THEN  PRINT  “  D«ri  ..»t  l  v«  Plot" 

1815  PRINT  L INC 2 > 

1316  IF  F1*g*0  Then  INPUT  “ENTER  MAXIMUM  ORDINATE  VALUE  Cmg>",C 

1817  0=0 

1813  IF  Fl*g*2  THEN  C*10O 

1818  IF  Fl*g*l  THEN  GOTO  1821 

1820  LINPUT  "DO  VOU  WANT  COMPUTER  CALCULATED  ORDINATE  T ICKS’" , Ax i t* 

1821  IF  F I »g*l  THEN  Axi*»*“Y” 

1825  A*FNM4*<XC» >,N> 

1330  B*FNMin<XC*'>,N) 

1835  IF  Axi*»tl,lJ<  "Y"  THEN  GOTO  I860 
1840  C«FNM»xCYC*1,Ni 
1850  D*FNMi nC Y<*> , N> 

I860  E*C  A-B  >.-■  10 
1870  F.CC-D)  10 

1371  IF  PI ot t «r *1  THEN  PLOTTER  IS  "GRAPHICS" 

1872  GRAPHICS 

1880  IF  Plott«r*t  THEN  LOCATE  28,90,20,90 
1835  IF  Plott«r*2  THEN  LOCATE  20.90,20,65 
1890  SCALE  B.A.D.C 
1900  AXES  2*E, 2*F ,  B.  D 
1910  MOVE  XC1>,Y<1) 

1920  FOP  1*1  TO  N 
1930  DRAW  X< I > , Y< I > 

1940  NEXT  I 

1950  MOVE  B,D-F 

1960  LORG  6  (LABEL  X  AXIS 

1970  FOR  I*B  TO  A  STEP  2*E 

1980  MOVE  I,D-.t«F 

1990  LABEL  USING  "DDD"JI 

20O0  NEXT  I 

2018  MOVE  B.D-F 

2020  DEC 

2830  LORG  8  t  LABEL  Y  AXIS 

2048  FOR  I*D  +  2*F  TO  C  STEP  2*F 
2050  MOVE  B- . 1«E,  I 

2055  IF  F 1 4g* 1  THEN  LABEL  USING  "D.DDE“lI 

2060  IF  F I 4g*0  THEN  LABEL  USING  "DDD. DD“ ; I 

2065  IF  F 1 4g»2  THEN  LABEL  USING  "DDD"; I 

2070  NEXT  I 

2030  LORG  6 

2O90  MOVE  8t3*E,D-.8«F 

2100  LABEL  G* 

2110  MOVE  Bt.4»£,D*5«F 

2115  LORG  4 

2120  L-DIR  90 

2  1  '8  LABEL  H* 

2131  LDIR  0 

2132  IF  P 1  ot  t  cr*2  THEN  OUTPUT  7,5;''IM" 

213:  LINPUT  "DO  VOU  WANT  GRAPHICS  OUTPUT  ON  THERMAL  PR  INTER" , Pri nt * 

2184  IF  Pr  int.»C  l,  1  ]*“•'“  THEN  dump  graphics 
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2133  II LEHR 
2136  £  IT  GROPmICS 

2139  .I-FUT  "DO  rOO  URNT  THE  PLOT  REPEATED"'” ,  ftp*  * 

2144  :p  Spt  *tl ,  1 1*'".-'  TriEN  GOTO  1301 
2149  i.SEXIT 
2130  3 . SEND 

2160  lEr  FNMi  <X<*.',N> 

2170  ••*  *X(l> 

2130  PIP  1*2  TO  H 
2190  IP  Mi  l  '  -Hit  THEM  I  > 

2200  ••£  T  I 

2210  RETURN  Mix. 

2220  FrEND 

2230  DEF  FNM i n<  7  *  *> ,  N  > 

2240  -in*Xtl' 

2230  FOR  1*2  TO  N 
2260  IF  XkI><Min  THEN  M-,  r*Xk  I ' 

2270  -E.T  I 
2230  RETURN  Min 
2290  P--END 
2300  3u8 

2390  R.jtEXIT 
2400  :U2END 
2500  t-EF  FMSeil'IX> 

2300  s.B  Spl  i r.4*  N,  Niro,  Xk * ' ,  Yi  4  . , Doinii r.‘  4 > , Func '  *> ,  Dir  i u< *>  #  Int ,  Ep» > 
2830  Ji3dti*iEpi  *0' 

2890  IP  Bid<lti«0  THEN  2990 

2920  'PINT  L  iNk  2  ,  "ERROR  IN  SOBPPOGRRM  Splini." 

2930  P  f  I  NT  ■Ep»*,,;Ep»,LlN<2  • 

2940  P4USE 

2930  GITO  2880 

2990  OPTION  BASE  1 

3OO0  tin  Si N> , G<  N-l > , Work •  N-l > 

3010  FOR  1*2  TO  N-l 
3020  Xi»XC) 

3030  Xt»l*Xk 1-1 » 

3040  Xi p 1*X< 1*1  • 

3030  Y  l  *Yk  I  ) 

3060  Yiit*Ytt-l> 

3070  pl*Yi 1 4 l  ) 

3080  '  *Xi«Xi»t 

309O  H*Xipl-Xi*l 

3100  Work. I)».5*X/H 

3110  T  *  ‘  <  V  i  p  1  -  V  i  )  -  •■  X  i  p  1  -X  l  1 - 1  Y  l  -  i  i  k.  1 1  .'X  <  -  H 

3120  S«I>«2*T 

3130  Gt I >*3*T 

3140  nEXT  I 

3130  :  1  >*SilO*0 

3170  **8-4*SQR<  3  > 

3180  v*0 

3190  FOP  1*2  TO  N-l 

3200  T»W*< -S'- 1  >  -Work  k  I  ■  I  - 1  >-<  .  3-Work  k  II  >*Sk  141  >  +G<  I )  > 

3210  H*R»S<T> 

3220  IF  HiU  THEN  U*H 

3230  Sk  I )  *S<  I  >  *T 

3240  -.EXT  I 

3250  I«  U>*Ep*  THEN  2180 
3260  for  1»i  TO  N-l 

3270  G' I >*<S' I»1 '-S< 1 1 • •  OU  l4i >-X‘  I  >  > 

3280  -EXT  I 

3290  Ie  Niro*0  THEN  3600 
3330  FOR  J* 1  TO  Nirj 
3340  lorrictor!  1*1 
3330  T*Doiii n>  J ' 

3360  IF  T  -iu  1  •  THEN  3430 

3370  PRINT  lIN>  2  i ,  "ERROR  IN  S'OBpFOGFRM  Selin*.'' 
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3380  PRINT  “ARGUMENT  OUT  OF  BOUNDS." 

9390  FRINT  TX  ••  1  X<  1  >,“'<■  H  •*"  j  X'  N/ ,  "Domai  ru  J;  “  =  Oiiuai  r.\ J >  .  L  1N>  2> 

3400  PAUSE 

3430  1-1*1 

3440  IF  I  N  THEN  3370 

3450  IF  T>XtI<  THEN  3430 

34o0  1*1-1 

3470  H*Dofc*intJj-y.<.I> 

3480  T*Doaai n< J >-X' 1+1 > 

3430  X*H*T 

3500  S*S<.  I  >*H*G'  1 1 

3510  2*1.6 

3520  U*2«(Sd)*SU*l)*S.' 

3530  W*<  Y<  1*1  >  -  Yt  l  >  >/<  X<  1  +  1  /  -Xt  I  }  '? 

3540  Func <. J)»W«H+Yt I >*X*U 

3550  D«r  i  u<  J>*U*(H*T>*U*2*X*G'.  I  > 

3560  NEXT  } 

360O  Int  *0 

3610  FOR  1*1  TO  N-l 

3620  H*.:<  1*1  >-X>  I' 

3630  Int*Int*.  5*H*<  Yi  I  >  +  Y<  I*  1  -1  '24*H  '3*t  S>  I  <*S<  I  *  1  >  ' 

3640  NEXT  I 
3650  SUBEND 

3800  SUB  K  i  n«t  i  c»(  A*.  i ,  Y<*>  .;<*>,  N,  Mz> 

3802  REM  X»TEMP,  Y-MASS,  Z*dM'dT 
3804  OPTION  BASE  1 

3809  DIM  Ff  tN>, T i nv<N> ,  XxtN> ,  Vy'NJ ,  D»CN> 

3814  INPUT  “ENTEP  RESIDUE  MASS  '  r»g>  “  ,  R*t  i  du« 

3320  Mr*Nz-R«t i duo  1  Total  mat*  r*»ct«d. 

3824  INPUT  “ENTER  POINT  NUMBERS  FOR  KINETICS  PLOTS” , P 1 , P2 
3829  Datapointa*P2-Pl*l 

3831  PEDIM  Ff  ■  Datapoi  ntj  >  ,T>nv<  Datapoi  ntt/ ,Xx'  Dat  apoi  n> »  >  ,  Yy  (Dat  apoi  nt  t  >  .  DhtDat 
apoi  mi) 

383*  FOP  1*1  TO  Datapontt 

3839  Xri  I  >»;:cPl  ♦!-!  1*273  !  Cowput*  array  and  Chang*  <o  K«luin. 

3342  V  y  <  I  ) *Mz  -Y r. Pl  +  I-1  >  '  Conput*  mast  r*act*d  array. 

3844  Ff ( I >»Yy< I >  'Mr  'Calculate  fraction  decomposed. 

3845  Ttnv< I >*1/Xxt I >  'Calculate  reciprocal  tthptrtturt. 

3846  NEXT  I 

3850  CALL  Spl  in*'  Dat  apo  i  nt  t ,  Dat  apo  i  nt  t ,  Xx<*  > ,  Ff  <  *  > ,  Xx(«>,Ff!.#>,Dm(»>,a2,l£-6.'lF 
it  cu rue  and  calculate  dCfracl/dT 
3355  MAT  Dm*AS$(Dm/ 

3863  CALL  PI ot z< Dat apoi nt s,X d ♦> , Ff <*> , D*<»> , A#> 1 Pr i nt  and  plot  dx.  dT  os  x 

3864  FOR  1*1  TO  Datapointt 

3869  IF  Dl»<  I  /  >.001  THEN  Dm'.  1  >*LOGC  1000*Dm(  I  >  > 

3374  NEXT  I 
3880  PRINT  PAGE 

3885  PRINT  A* 

3886  PRINT  LIN<2> 

3890  PRINT  "Lin**r  r*gr»t*ion  analysis  of  TGA  data  between  Temperatures" 

3980  PRINT  USING  3910; Xx< 1 > , Xx< Dat apoi nt t> 

3910  IMAGE  DDDB, ”  and  *,DODD|“  d*gr**«  Kelvin. * 

3920  PRINT 

3930  CALL  Regress.  Tt nv(*>,Dm< *>,Da> apoi nts, A, B,C),C2,C3, C4,C3,C6,C?,C8, C9, CO,  Va 

r  > 

3940  PRINT  USING  3950; ABS< B*1 . 9872) , Van 

3950  IMAGE  “CALCULATED  ACTIVATION  ENERGY  IS  ” , MD. DDE , *  Cal.xaol*  *,«  ♦/-  “,D.DP 

3951  PRINT 

41 10  PRINT  PAGE 
4120  PRINT  LINC3) 

4130  PRINT  At 
4140  PRINT 

4J50  PRINT  “KINETICS  DATA  PLOT" 

4160  PRINT  LIN<2> 

4170  x»a;  *FNMax < T ■  nv( *  • ,  D a’ apo i nt  »  > 
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4190  .  :»t  n*FNMi  n<  Ti  n>  \*) ,  t  »t  *po  1  r.t  *  > 

4190  Vtii*x*FNM».~  <  Bm<  *  .> ,  B»t  apo  i  nt  *  > 

4200  Van  r>«FNMi  ncB«<  *  >  ,  B*t  tpo  I  nt  i  > 

4210  E*vX»»x-Xmi  n.'-- 10 
4220  F-i'YM»x-YMin> '10 

4222  INPUT  "ENTER  PLOTTER  COBE:  1-CRT;  2*9872A" . Cod« 

4224  IP  Cod«-2  THEN  PLOTTER  IS  7,5,“98?2A" 

4230  IP  Cod**l  THEN  PLOTTER  IS  “GRAPHICS" 

4240  GRAPHICS 

4250  LOCATE  10,90,15,75 

4260  SCALE  Xm i n, X»*x , Van  n-F, Ym*x*F 

4270  AXES  2*E,2*F,XMin,Y»in-F 

4280  FOR  I»1  TO  B*t*point» 

4290  MOVE  T i nv< I > , Bn< I ) 

4300  LABEL  “+“ 

4310  NEXT  I 

4320  FOR  J»1  TO  Batapoint* 

4330  B»i J /*B*T i nu< J ) +A 
4340  NEXT  J 

4350  MOVE  Tinu(P,l*(P 
4360  FOR  I»1  TO  Batapoint* 

4370  BRAN  Tir»v<I>,BM<I> 

4380  NEXT  I 

4398  HOVE  Xialn,  Yiain-F 

4400  LORG  6 

4410  FOR  I*X«,n  TO  Xmax  STEP  2*E 

4420  MOVE  J,Y«in-t.l»F 

4430  LABEL  USING  "B. BOB"; 1*1008 

4440  NEXT  I 

4450  MOVE  Xa.in,Ya>ir>-F 

4460  LORG  8 

4470  FOR  I«Y»irv*F  TO  Yiaax+F  STEP  2*F 

4480  MOVE  X«in-.04»E, I 

4430  LABEL  USING  "B.BDB“;t 

4500  NEXT  I 

4510  LORG  6 

4520  MOVE  X»in*5*E,Y«iin-2*F 
4530  LABEL  “1/T<K>  x  1000" 

4540  MOVE  XMir>*.4*E.YMin*5*F 
4550  LDIP  PI'2 

4560  LABEL  “In  >dM'dT>  x  1000' 

4570  PAUSE 

4574  LINPUT  “BO  YOU  WANT  TO  REPEAT  THE  PL0T7“,Rpt* 

4579  IF  RptBCl, 1]«“Y“  THEN  GOTO  4222 
4534  IF  Cod««l  THEN  BUMP  GRAPHICS 
4594  EXIT  GRAPHICS 
4608  GCLEAR 

4604  LINPUT  "BO  YOU  WANT  TO  REPEAT  USING  BIFFERENT  PO INTS?" , Agat n* 

4608  IF  Ag»tn»tl, 11*“Y“  THEN  GOTO  3824 
4610  SUBEXIT 
4620  SU1ENB 

4700  SUB  R«gr**»<X<*)fY<*>,N,A.B,R«gJ»,P4***,Tot*l»», R*g»» , R»*»* , F ,Bfr«g,Dfr*»,B 
f  tot , Abort , B«v> 


4710 

■  •*♦  mobel:  y«a*b*x  ««* 

4720 

ON  ERROR  GOTO  Boab 

4730 

Abort  >0 

4740 

Yl«Xl*Z»X2»Y2-0 

4750 

FOR  I«1  TO  N 

4760 

X1*X1+X< 1 > 

•  Sup 

of 

X't 

4770 

Yl-YWU  ) 

!  Sum 

of 

Y'» 

4780 

X2»X2*X<:I>*X<I> 

!  Sum 

of 

X  *quar«* 

4790 

Y2-Y2*Y< 1 > *Y  < I ) 

i  Sum 

of 

Y  »qu»r»i 

4800 

Z*Z*X< I )»Y<1 > 

•  Sum 

of 

XY'» 

4810 

NEXT  I 

4820 

X1«XI  N 

>  M car 

1  X 

4830 

Y 1 «Y 1  N 

!  M*»n 

1  Y 
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•*340  $*•  i-N*Xl*Yl ) 

4350  A*Y1-B»X1 

4860  Tot il ss*Y2-N*Y14Y1  !  Totil  Sum  of  Squirts 

4S70  Reg**miZ-NeXl»Yl >*2  <X2-N*X1«X1 > i  Regression  Sum  of  Squires 
4330  Resss«Tot il ss-Regss  !  Residuil  Sum  of  Squires 

4393  R«gms«Regss  1  Regression  Hein  Squires 

4983  Res«is«Re*ss-' (M-8 '•  !  Residuil  Hein  Squires 

4904  De'-,»SQR<Resms /  1  Mein  Squire  Devi  it  ion 

4310  3FF  ERROR 
4928  DEFAULT  ON 

4930  F»Regms  Resms  !  F  Rit t o 

4940  DEFAULT  OFF 

4950  Dfreg»l  !  Degrees  of  Freedom 

4963  Df re j  «N-2 

4970  Df tot*Dfreg+Dfres 

4930  SUBEXIT 

4990  Bomb:  Abort *1 

5000  SUBEND 

5330  SUB  Exponent  i  il  t  X<*) ,  '/(.*  > ,  N,  A,  B,  Regss,  eesss,  Tot  it  it ,  Regms,  Re  sms  ,  F,  Dfreg,  Df  r 
es , If  t  ot , Abort ) 

5340  I  model:  Y»A»EXP<B«X>  *** 

5330  ON  ERROR  GOTO  Bomb 
5360  Abort »0 
5370  X1«X2*Z*Y1»Y2»0 
5380  FOR  I«1  TO  N 
5390  Logy”L0G<  Y< I ) ) 

5460  X1»X1*X<I>  !  Sum  of  X's 

5418  X2“X2+X< l >»X< l )  1  Sum  of  X  squires 

5420  Yl-Yl+Logy  <  Sum  of  Y's 

5430  Y2»Y2+LogyitLogy  !  Sum  of  V  squires 

5440  Z«Z  +  X<  I  )«Logv  !  Sum  of  XY'S 

5450  NEXT  I 

5460  Xl»Xt/N  ‘  Mein  X 

5470  Yl-Yl/N  1  Hein  Y 

5480  8*1'  Z-N»XI*Y1  >  '<X2-N*X1*XI > 

5490  A»EXP<Y1-B*X1) 

5500  Totil  ss”Y2-N*Yl»Yl,  I  Totil  Sum  of  Squires 

5510  Regs4«<Z“N*;<l»Yl)'2'<X2-N*Xl»Xl)!  Regression  Sum  of  Squires 
5520  Resss»Toii’ss-Regss  1  Residuil  Sum  of  Squires 

5530  Regms*Regss  !  Regression  Mein  Squires 

5548  Resms«Resssr<N-2)  1  Residuil  Mein  Squires 

5550  OFF  ERROR 
5560  DEFAULT  ON 

5570  F«RegmsrResms  1  F  Rltto 

5580  DEFAULT  OFF 

5590  Dfreg»l  !  Degrees  of  Freedom 

5600  Dfres«N-2 

5610  Df tot »Df reg*Dfres 

5620  SUBEXIT 

5630  Bomb:  Abort >1 

5640  SUBEND 

5650  SUB  RIoti<N,T<*>,X<*>,Y<»>,A»> 

5660  OPTION  BASE  1 
5670  PRINT  PAGE 

5680  PRINT  LINO) 

5690  PRINT  A* 

5695  L INPUT  “DO  YOU  WANT  THE  DATA  PRINTED?", Print* 

5696  IF  PrinttCl,lI»“Y“  THEN  COSUB  Print 

5700  PRINT  "Kinetics  plot  of  dx/dT  us.  x.  *■  friction  decomposed" 

5710  PRINT  LINO) 

5720  C«FNMix<Y<  *  >,N> 

5730  D*FNHt  n ' Y  ;e  > , H) 

5740  E». 1 

5750  Fe<  C-D  10 

5754  INPUT  "ENTER  PLOTTER  CODE:  1-CRT;  2-9S72A  ".Code 
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5760  IF  Code*l  THEN  PLOTTER  IS  "GRAPHICS" 

5764  IF  Code*2  THEN  PLOTTER  IS  7,5,"987ZA” 

5770  GRAPHICS 
5780  LOCATE  25,90, 15.75 
5780  SCALE  O.l.P.C 
5800  PTES  . 2, 2*F ,  0,  D 
5810  F OP  1*1  TO  N 
5S20  MOVE  X-n.Y'l) 

5330  LABEL  “♦* 

5840  NEXT  I 
5350  LORG  6 

5360  FOR  1*0  TO  t  STEP  .2 
5870  HOVE  I,D-.1*F 
5880  LABEL  USING  “2. DOD“ ; I 
5890  NEXT  I 
5900  LORG  8 

5918  FOR  I*D*2*F  TO  C  STEP  2*F 

5920  MOVE  -.01,1 

5938  LABEL  USING  “D. DDE" I  I 

5940  NEXT  I 

5950  LORG  6 

5960  MOVE  .S,D-.3*F 

5970  LABEL  “Fraction  Dec ompostd" 

5990  LD1R  PI/2 

6000  MOVE  .823,0+5*F 

6018  LABEL  "Rat*  <»g<-deg.>" 

6020  LOIR  0 
6030  PAUSE 

6040  IF  Code*l  THEN  DUMP  GRAPHICS 
6850  GCLEAR 

6054  L INPUT  -DO  YOU  UAH'  TO  REPEAT  THE  PL0T7",Rpt* 

6055  IF  RpiBtl,  11*“V  THEN  GOTO  5700 
6060  EXIT  GRAPHICS 

6070  SUBEXIT 

6000  Print:  PRINTER  IS  0 

6098  PRINT  "TeaperatureOO  Fraction  decoapoaed  Pate  Ideg^-1)* 

6100  PRINT 

6110  FOR  1*1  TO  N 

>120  PRINT  USING  6150;T<I>,X<1>,YU> 

6130  NEXT  I 

6150  IMAGE  DDDOXXXXXXXXXXXX,Z.DDDDXXXXXXXXXXXXXXXXXX,D.DDDE 

6160  PRINT  PAGE 
6170  PRINT  LINO) 

6180  PRINT  A« 

6190  RETURN 
6200  SUBEND 
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Appendix  B 

Program  for  Acquisition  of  Isothemal  DSC  Data 
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FRGSF  FOR  ■=!  ■  J I S  ITION  An?  ANALYSIS  OF  ISOTHERMFiL  DSC  DAT  ft 


10  OF*! OH  BASE  1 

20  j:-  V  1000' .  Turn  1000  ,:3<1000  )  ,V3>.  1000s 

30  printer-  is  is 
40  PR  I'-T  LIN'  4 j 

50  PF  InT  "lift  aquijition  progr  am  for  th*  Scanning  Bi  t  i  al  " 

SC  PSI*IT  "CtloriM’tr.  U»*  scanning  digital  1  i  n«  " 

70  PF  I  NT  "0.  Ordina**  lxit  i*  calibrated  for-  10  Millivolts" 

S0  P'lNT  “full  *cal*.“ 

90  PF  :*1T  "Da*a  i  i  tal*--.  »i  th#  approximate  rai  e  of  or,*  point” 

100  —  INT  "P«f  4  *«co-'ii.  Edit  Mn«  313  to  chang*  rat*" 

110  F  4  I  NT 

1D0  “  *  t  N  T  "Follow  *  h«  CPT  pr  jt.pt-  s  •  Pfttt  CONT  to  *nt«r  data" 

110  ff  INT  “into  th*  c  oiMput  *r  or  to  haw*  a  plot  on  th*  CPT  put" 

140  PRINT  "onto  th*  thermal  printer.  Terminate  th#  progr  or*  by  preiiing1' 

145  rt INT  “STOP." 

146  PRINT 

150  =  =  INT  "EDIT  THE  PROGRAM  AND  PRESS  CONT" 

155  CAUSE 

160  =  f  I  NT  PAGE 

170  PSINTER  IS  0 

ISO  I IM  AfC60],BfC602.C*C803 

ISO  -INPUT  •  ENTER  TEST  IDENTIFIER,  LINE  1“,A* 

200  -INPUT  "ENTER  TEST  IDENTIFIER,  LINE  2“,B* 

210  I-.FUT  "  ENTER  TEMPERATURE  OF  THE  RUN",  Ti  nit 
225  INPUT  "ENTER  MASS  C'F  SAMPLE  (.mg)  "  ,  Samp' *maf  » 

220  INPUT  "ENTER  SENSITIVITY  >  meal  r»#c  >  * ,  S*n»i  t  i  vi  ty 
240  Pc INT  LIN' 3 i 
2S0  POINT  A* 

260  POINT  B* 

2*0  PRINT  USING  2901 T i ni t 

2’5  PRINT  USING  285; Sampl email, Senai t  ivi  ty 
283  I ''AGE  "Temper  at ur*  il  “»DDD,"  K  “ 

285  IMAGE  "Mail  of  lample  ij  '.DD.DDD,"  ng.  S«n»i'i"Uv  it  ".DD,"  mcal  tec” 
290  PRINT  LIN. 2) 

300  OEM  Th*  following  tin«i  configure  th*  scanning  digits'  voltmeter  and  tin* 
clod  and  s«t  v«y  int*rupts. 

310  OUTPUT  701;  "R  h2  0"  'Configur*  DVM  for  inurup'  'ransf*". 

320  l!N  ERROR  COSUB  Etcap*  1  Enab  '  *  run  »rrors  to  b*  corntua. 

3S0  KEY  »7  GOSUB  End'Op*rator  t*ri»i nation  of  data  aquuittcn, 

335  3N  KEY  »8  GOSUB  Exit 

340  3N  INT  #7  GOSUB  R*ad  '  Br  anchi  ng  routin*  for  data  tranif  «r. 

350  OUTPUT  ?; "ArUl*I2“ IConf igur*  clock. 

360  I0NTR0L  MASK  7;  128'Look  for  int«rupt.  flag. 

370  'PINT  L I N< 2 > 

375  PPINT  "Ti  a*'  i*c  ■  Ordi  nat*  Pol  nil1' 

380  FP INT 
3S5  1*1 

390  DISP  "Pr*ss  CONT  to  b*gin  data  aquuition" 

400  »AuSE 

410  OUTPUT  9; "U1G" 'Start  iht  clock 

420  IISP  "PRESS  t7  TO  TERMINATE  DATA  GATHERING" 

430  CARD  ENABLE  7!Allow  DVM  int*rupl*  to  tak*  plac*. 

440  if  i<ieee  then  goto  440 

450  DISABLE 

460  PRINTER  IS  16 

470  PRINT  "ALLOCATED  CORE  STORAGE  IS  FULL" 

480  PRINTER  IS  0 


NWC  TP  6337 


4  30  OOSUB  End 

500  E  > t :  OUTPUT  701 ; 

505  3TO» 

510  P«»C!  '  INPUT  IATA  FPOM  DSC 

520  entef  roi ;  •  •  I  -• 

52?  V  I  »V'  I  *1000+20  ‘Co-'./trt  Millivolt*  to  a  3  to  10  i.'llt. 

530  OUTPUT  * ; “ U 1 V " 

540  ENTER  8;  Ti 

550  T  i  n,*  •  I  **T  i  me  •  I  >.  1008  'Convert  ui  1  I  iitcoridi  tc  second*. 

560  PRINT  USING  370:  Ti  me  •  I  ■>,  Y<  I  <,  I 

570  IMAGE  DDDDXXXXNXXXXVXX.MDD.DDXXXXXXXX.DDDD 

530  I*I»1 

630  OUTPUT  70i;"0"  'Trigger  Digital  Voltmeter 
660  CAPD  ENABLE  7  'Re-Enabl*  I  nt  erupts 
670  RETURN 

630  End:  '  Pro:***  Lit* 

685  INPUT  "ENTER  POINT  NUMBER  OF  ST  ABLE  ORDINATE" ,  R 
687  FOR  1*1  TO  R-l 

630  -I  ■«-'<P>  >El  KMnv  i  initial  transient 
605  next  i 

700  hpis*I-2  'Eliminate  jxirintoui  dlt*  point* 

710  PEPIN  V1'  Npt  *  i ,  Ti me <Npt * /  'Cond«n**  the  arrive. 

720  Ml ni mum»FNMi n' Y- « > , Npt * >  'Ca’culate  I  outs  i  voltage  output  Fro*  PSC. 

730  MAT  Y»Y-> Mi  hi  »u»,'  !Se*  jtro  point  is  lowest  output  voltage 

’70  CALL  Gr  aph  >  A*,  B*,  T  i  me  <  *  >  ,  Y<  ♦  ' ,  Npt  s')  I  Plot  the  dl*  *. 

780  LIN01.IT  “DO  YOU  WANT  THE  DATA  STORED'’" ,  C* 

730  IF  C»Cl,ll«"Y“  THEM  CAlL  St  or*  ■:  A* ,  B* ,  T  i  me  t  ♦  > ,  YC*  > ,  Npt  *  ) 

800  LINPUT  "DC  YOU  WANT  TO  EXPAND  THE  PL0T",D* 

310  IF  D*Cl,13»"Y”  THEN  CALL  E  p»nd< A*, B*.  Ti »*<♦ > ,  Yt *3 ) 

820  LINPUT  DO  YOU  WANT  TO  MATE  A  LINEAR  BASELINE  CORRECT  ION? ", E* 

830  IF  Eltl,l>’V  THEN  CALL  Trin*for*,<A*,B*,Ti**<*>.Y<*i,Npt»> 

335  FRINTEP  IS  16 

836  PRINT  PACE 

337  PRINT  "P«al  integrator  using  in  arbitrar"  baseline  correction" 

833  PRINT 

838  PRINT  "Data  point  number e  correspond  to  origional  cat  a  list" 

340  PRINT  LIHiI, 

341  PRINT  “The  baseline  ie  calculated  using  a  jtught  line  from  the  firs*  to 
the  last  data  point" 

842  PRINT  LINc3> 

843  PRINT  "Press  STOP  if  you  wish  to  terminal*  th«  program.  Press  CONT  if 
you  wish  to  continue." 

844  PAUSE 

845  PRINT  PAGE 

84?  INPUT  'ENTEP  BEGINNING  AND  ENDING  POINTS  fop  NUMERICAL  INTEGRATION". Begin, 
End 

850  Opt s*End-B*gi  n*l 

852  PRINT  USING  854;0pts 

.854  IMAGE  "ENTEP  1  AND  ",DDD,"  FOR  INTEGRATION  LIMITS" 

856  REDIM  X3tCpts  ■,Y3<0pt») 

836  FOR  L»1  TO  Opt* 

860  Y3' L.'«Y' L*Begi  n-1  > 

862  X3'  L  ■»Time',L+legir.-l  > 

864  NEXT  L 

866  SI cp*3* < >  3 ' Opt  s  > -Y3 c 1 1 > • < X3< Opt s  > -X3 • I > ) 

867  B*Y3< 1 > 

868  FOP  L*1  TO  Opt* 

870  Y3<L>»Y3tL>-Slope3*tX3tL)-X3<l'»-B 

372  NEXT  L 

874  C»*I*4” :  Expanded  plot  peak  integration" 

376  CALL  Graphic  A*  ,CP.X3<*>,Y3'.  *)■  Opt  a.  Area: 

878  LINPUT  "DO  YOU  WANT  TQ  TERMINATE  THE  PPOCRAM 7" • St  op* 

880  IF  StopfC 1 .  1 v*  THEN  GOSUB  Exit 

881  CALL  1 1  net  i  c  s  :.3<*  > ,  Y3t  ♦) ,  Oot  s,  Ar e  i’ 

884  GOTO  300 

833  Escape  1  Correct  any  run  errors. 
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s?o  -a  inter  is  is 

200  -=:  *T  ERRmi;."  GENERATED.  EDIT  PROGRAM  AM'  PRESS  coht.- 

*10  rause 

*20  FR INTER  15  0 
930  :ET..RN 

.000  :  JS  St  ore'  ft*.  B*.  *  < ,  I  .  t  • ,  M  > 

1010  r  =  TJOM  BOSE  1  1  SuS/ou*  in*  t:r  I'cr  i*;  3*'  a  fro*.  '  ~V  plotter  input. 

1011  SHORT  UUO.V'.N' 

1012  ”AT  u»>: 

101 3  HAT  u«Y 

1020  -INPUT  'Enter  MOSS  STORAGE  MEDIUM ',r# 

10  30  -InF-jT  "ENTER  i  CHARACTER  FILE  NAME",GI 
1040  ■::>** 

1050  CREATE  C».N  9 
lOCO  ASSIGN  #2  TO  Cl 
10T0  RR I  NT  »2;AI,B*.N 

ioso  hat  print  »2; v 

1090  F  f I  NT  LIN- ?  • 

loss  -rint  using  i ioo: t= :mi  gi 

HOC  : "DATA  IS  STORES  IN  FILE  "  .  R  ,  “  " 

1110  ASSIGN  *2  TO  » 

1120  SUIEYIT 
1130  SJBEflD 

1200  SOS  Graph<A*.B#,X'« -.Vt*>.N> 

1210  RRIS*  PAGE  ’$L9.'0ut  in*  for  plotting  ga*a  from  an  X-Y  plotter 

1220  PRINT  LIN/3i  ’tiling  the  internal  ther**)  printer. 

1220  OPTION  BASE  1 
1240  PRINT  Al 
1250  RRtHT  B* 

1 200  RRINT  HN.'3> 

1220  2IM  DIC19] 

122?  C*»"TIME  •  n:  •“ 

1220  t'l«“Ri  1  4*  '  '*  M  !  1  1  l  vol  t  i  " 

I  200  F*FNM»'  <Y‘  ♦  .  ,  N,  1  • 

1 290  S-  *'i*.  N  > -X 1 1  * 

1300  A«F 

1310  -INPUT  'to  YCU  WANT  n  MAYU’Uh  Y  USES  FOR  THE  GRAF m  INSTEAD  OF  THE  CALCULAT 
ED  MA.  IMOM*,J* 

1  320  IF  J»tl.ll«-Y*  THEN  INPUT  “ENTER  Yi*a*”,H 
1325  IF  TICI. 1]«"Y"  THEN  A»H 

1  330  v INPUT  "ENTER  PLOTTER  CODE:  CRT-fl.  9322A*B " ,  P  i I 
1S40  IF  Flotici.  13*"B"  then  '-OTTER  IS  2.5,  "RSTCA’' 

1350  CF  P'otICl,  11»"A“  THEN  PLOTTER  IS  ''GRAPHICS" 

13E0  2RARHICS 

li'O  -OC ATE  IS, 33,10.2? 

1 3SO  SCALE  :!<  1  i.X'N  ■  .O.A 
1330  A  YES  B-4.A  10, X<  1  >.0 
1400  MOVE  X  'll ,  Y  ••  1  ■ 

1410  FOP  I»1  TO  N 
1420  DRAW  I  >  ,  Y <  I  » 

1430  NEXT  1 
1440  PENUR 
1480  -ORG  i 

1490  FOR  !*-;<»  •  TO  :;.N'  STEP  8/4 

1?0O  N'jYE  I. -A- lOO 

1? 1 0  uAbel  using  i?3o:  1 

1 ?20  NEYT  I 
1 ?30  IMAGE  DDD 
1  ?40  LORO  9 

1?50  FOR  I. A  10  TO  A  STEP  A/10 

l?eO  MOVE  -B-lOOft  i  ',1 

1520  -ABEL  USING  1590:1 

1590  NEXT  I 

1530  IMAGE  I'D.  D 

1000  -ORG  i 

1C10  'JOVE  B  2-Y  t'.-A  '.? 
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1620  L  A E  L  CS 
1630  '-DIP  PI'2 
1640  LOPC  4 

1650  MOVE  M  ■,d  '2 

lose  label  d* 

1670  w INPUT  -BO  VOU  WANT  GRAPHICS  PRINTOUT  VIA  THE  THERMAL  PRINTER?" , L* 

1630  IF  Lit ! . 1 J*"V"  THEN  DUMP  GRAPHICS 
16*0  E:  IT  GRAPHICS 
1700  GCLEAR 

i ~os  if  Pio*iti,  n="B”  then  output  res;  “Spc- in- 

1710  SUBEXIT 

iree  subend 

1730  DEF  FNMaxv  A1  4  1 ,  N,  St  ar  t  " 

1740  M*,.=A>.  St  art  > 

1750  FOP  I  •Start  TO  II 

1760  IF  A'  I3  Max  THEN  M»x=Av  I  > 

1770  NEXT  I 
1730  RETURN  Mis 
17*0  FMEnt 

3000  SUB  £  'pandICI|t#,IX»>,V(«)' 

3010  OPTION  BASE  1 

3030  DIM  3>  1000', T<  1OO0>  ,  R<  1 000' ,  «<  1000.' 

3035  DIM  I*C30J,P#C3Q3 

3030  INPUT  "ENTER  LOWEST  AND  HIGHEST  DATA  POINTS  FOR  PLOT  EXPANSION  •.  PI .  F2>  “ .  PI 
,P2 

3040  M«ABS iPi-Pl >♦! 

3050  RED  I M  2<N>,T<N> 

3060  FOR  1*1  TO  N 
3070  ;>  i  ■•x',i*pi-i> 

30SO  T.  I  ..Vtl-Pl-n 
30*0  NEXT  I 

3100  I  f  «D#C.  “  1  ExparOat)  Plot” 

3110  ‘.INPUT  “IS  THIS  A  BASELINE  CORPECTED  PLOT  ,  M* 

3130  IF  M*tl,13-'x“Y”  THEN  GOTO  3175 

3 1  SO  IF  n*C  1*1  3*“ Y"  THEN  CALL  PnnttlCI,  IP,  3c*  .  T *  *  1  ,  N 
31*0  IF  M*;i,n«'V“  THEN  CALL  Cr  »ph2<C$,  I»,  Z<  *  J .  Tt  ♦  ' .  N,  Ar«  a) 

3141  PRINTER  IS  16 

314-  PRINT  "DO  YOU  WANT  TO  INTEGRATE  THIS  EXPANDED  PLOT  USING  ANOTHER  BASELINE 
CORRECTION'" 

3144  LINPUT  N* 

2145  PRINT  PAGE 

3146  IF  N«C!,n'  ,“Y“  THEN  GOTO  2170 

214-  INPUT  "ENTER  THE  DATA  POINT  NUMBERS  •  taw  list'  FOP  ORDINATE  POINTS  ",P3,P4 

214;  Mx *F 4-F  3* 1 

214-  FEDIM  ft'Mxl,N(Hv) 

2150  PPInT  USING  2153;Mx 

2152  IMAGE  "USE  1  AND  ” , DDD, ’  FOR  INTEGRATION  LIMITS" 

2153  FOP  1*1  TO  Mx 

2154  p/f »«2<K*F3-l> 

2155  W<2  .«t,k»P3-1' 

2156  NEXT  t 

2157  Rf«'E*P  irdad  Plot:  Arbitrary  baialina  correction" 

2153  S I  *p«2*'  U<  M  i  >-Ht  I  >  '•  <  RviNx  j-R<  l  )  ) 

215*  3b*U'  1  • 

2162  FOR  '•!  TO  llx 

216:  W't  '«U't  --S I  epa 2*t  P1-  K >-Rt  l  J  )-Bb 

2164  NEXT  2 

2165  CALL  Gr»pK2',C*.R*,R<»>.W(»),My» 

2170  IF  im,lI*"Y'  THEN  SUBEXIT 
2175  call  CrapnCC*.  It, 2t*),T<»),N) 

2  I  SO  S1JBEX  I  T 
2135  SUBEND 

2160  SUB  Tr  am  f  or  n>  F#t  G» ,  Xc»> ,  V  •'  *  ■ ,  N) 

2200  CST ! CM  BASE  1  1  Linear  baseline  co meet -on  For  2  -tenable  plot*. 

2210  DIM  Vr.«*i  N' 

22I5  DIM  H*  C  SO  1 
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.  7  Wi  •  3  .♦  •  9*Y,;0  >  5  !  A,"ir4.j»  *h  »  «;«•»*  4r.d  ’»i’  5  CC'fJ  »r9 

t  4k  «  'riiikl  dip  into  : :  t :? ’d -5 r’  «l  '  on 

2250  w»«.  y.  s  •*’•••  H-l  »*  •  •  N-2  •  .4-3  *>  P-4  '  5  1  'or  »!cp«  :  4’ •:  :j  1  »»  i  e«. 

2240  5  I  G  p  €  =  -  V  -  5  •  *  N  •  -  1 

2250  p op  i«i  to  n 

2160  .  r.*w  I  ■*  c  >  I  '-S '  op*  *•  I 

21"  ne:  t  1 

p*«g*t  " :  l  i  '•«  4-’  Suii'm  Ccp-s:'  i  on' 

22  ?0  C*i _ jt  4i:'G  fl.HI,  -  >,  N.  Ar«  4  > 

1550  lINF'.  t  'io  .00  hhNT  TWE  CLCT  EXPANDED’” .  iU 

133?  ;r  M*: ; , i :«"Y"  then  call  e  p4r.d' e*.h», ^ ■  -» • 

1340  -tJEE'.IT 
1350  SuBEND 

1400  DEE  ♦  .N.- 

1410  M 4 -  =  ABS>  A'  1  ■  > 

1410  F OP  1-1  TO  N 

1430  :E  *BS>A'  I  >  '  Mk  4Ci  THEM  «4- 4bi*0Bi'  0'  I  >  • 

1440  NEXT  I 

144  1  RETURN  fill  4b  i 

1441  RNEnD 

2450  2EF  run  1 IV  &■  *  • ,  N  • 

1460  Mi n*S< 1 1 

1470  EOF  1*2  TO  N 

2430  IF  B<I.'  Mir,  THEM  Mm*BI.' 

2490  NEXT  I 
25O0  RETURN  a  in 
1510  FliEND 

1600  SUB  i>kPPl  0*,I*.X' ■  *  1 .  N,  in  €  .;r4l  I 
160?  PRINTER  U  0 

1610  FEINT  EhGE  '3'jtro..'  m.  r'-jr  pi-jtting  d4*.4  t'oi.  4-i  •  ->  pl:*t4r 

1620  FPINT  LINO'  'v.4--r„3  -he  ir-**rn*l  •b*rn4l  printer  . 

1430  OPTION  BASE  1 
2640  PE  I  :iT  A* 

1650  PPINT  B* 

1660  PRINT  LIN'  5 • 

1670  DIM  I'*C  193, 0.  •  1000  ••  Vy-  1000  • 

267:  REDIM  U.\'  N0,V-."  N> 

1675  C**"TtME  'ha-" 

167S  Ii»*  P«  '  41  i  >•«  Mi  I  I  i '<c- 1  >  i" 

1650  E  =  F,.M4-4bi'  '■  *>.N> 

26  ?0  2*/  ■.  ■'  l  1  -  X  1  > 

2700  A*F 

2”10  LINE1.*  "DO  YOU  WANT  9  MAX !  MU''  USED  FOR  THE  GF*Pn  INSTEAD  OF  THE  CAuCUlAT 
ed  nfirinun* , J* 

2710  IF  J*Cl,  l  ]  *  "  >  "  Then  inFl-T  "Enter  .r«4,<",H 
2740  IF  .'*11,13=  V"  THEN  A  =  H 

2745  LINPOT  "ENTEF  PLOTTER  CODE:  CRT=A,  9E72fl  =  E- " .  P 1  at  s 
2750  IF  P1C'**tl.  U»“A”  THEN  PLOTTER  IS  “GRAPH  I  OS " 

1755  IF  PlotfU,  1  J*“B"  THEN  PLOTTER  IS  7,5, "55720" 

2760  GRAPHICS. 

2770  LOCATE  17,92.10.75 
2730  SCAlE  X  •  1  > ,  X  M.-A.A 
2790  AXES  B  4,  ft  •  10,  '•.<  I  ■ .  0 

1300  MOVE  :■  1  1  • 

2310  POP  1=1  TO  N 

2320  DRAW  ■  I  I  ■ 

1350  NE/T  I 
1340  PENUP 
1330  LORO  6 

2390  FOR  I*Kili  TO  SMI'  STEP  B  4 

2900  MO'-'E  1,-0  100 

1910  LABEL  USING  1930;i 

2920  NEXT  I 

29  30  I  MACE  DDD 

1940  wORC  3 

1950  F C P  I*-4  TO  A  STEF  A  !0 
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2960  MOVE  -B'  lOO*X<  1  • ,  I 
2970  LABEL  USING  2990:  I 
2920  NEXT  I 
2990  IMAGE  MDD.D 
3000  LORG  6 

301O  MOVE  B  2  +  X*l  *.-A  15 
3020  LABEL  C* 

30 SO  lDIR  PI  2 
3040  LOFO  4 

3O50  MOVE  -E  ?♦::<  1  ',0 '2 
3060  LABEL  El 

3062  LINPuT  "DO  YOU  WANT  TO  INTEGRATE  PART  OF  THE  CURVE", 0* 

3065  if  o»n,iJ,.  "r  then  goto  3035 

3066  INPUT  "ENTER  BEGINNING  AND  ENDING  DATA  POINT  NUMBERS*: FROM  NEW  SET  IF  EXPAN 
BED >“,St art, Finish 

3062  REDIM  Uxt F i ni sh-St  art  + 1  * ,  V; * i  F i  n i  s h-St  art  + 1  >  t 

3063  3*1 

3063  FOP  I  *  S  t  ar  t  TO  Finish 

3070  U  \  *.  J  >*X<  I  > 

3071  Vy  <  J  > *Y <  !  > 

3072  MOVE  I  >,'<■<  I’* 

307  3  DRAW  '  I  >  ,  0 

3074  PENOP 

3075  J*J»l 

3076  NEXT  I 

3077  Epsilon*lE-6 

3078  CALL  Intunt  J-l,Ep*i  1on,U:.<'«?,  Vyt4>,  Int«gr*l  3 
387?  PRINT  USING  308O;Uxt l>,Ux< J-l ), Integral 

3030  IMAGE  “THE  AREA  REPRESENTED  BY  THE  SHADED  REGION  BETWEEN  DOMAIN  VALUES  " , 
DDDD.DD,"  AND  ", DDDD. DD. “  IS  " , MD . DDDE 

3031  PRINT  LIN<23 

3035  LINPUT  "DO  YOU  WANT  GRAPHICS  PRINTOUT  VIA  THE  THERMAL  PRINTER7" , L* 

3090  IF  L»Ct,13*"Y“  THEN  DUMP  GRAPHICS 
3095  EXIT  GRAPHICS 
31O0  GCLEAR 

3105  IF  Pl6tm.lJ«"B"  THEM  OUTPUT  7,5;  "SPO-  IN" 

3110  SUBEXIT 
3120  SUBEND 

3130  SUB  Pnntz<ft»,B»fX<*>,V*  **,N> 

3140  OPTION  BASE  1 

3145  LINPUT  "DO  YOU  WANT  THIS  PRINTED  ON  THE  Ti PEWP I  TER?" . Rr i nt«r* 

3146  IP  Pri  nt«r  *[  1 ,  U*"Y”  THEN  PRINTER  IS  7,4 

3147  GOTO  3160 
3150  PRINT  RAGE 

3160  PRINT  "DATA  POINT  RENUMBERING  FOR  BASELINE  CORRECTED  .  EXPANDED  PLOT" 

3170  PRINT  A* 

31*0  PRINT  21 
3190  PRINT 

3200  PRINT  "DOMAIN  VALUE  RANGE  VALUE  NEW  DATA  FQINT  NUMBER” 

3210  PRINT 

3220  FOR  1*1  TO  N 

3230  PRINT  USING  3250; X< I  • , Y <■ 1' ,  I 
3240  NEXT  I 

3250  IMAGE  DDDD .  DDXXXXXXl.XXXXX,  DD.  BDXXXXXXXXXXX,  DDD 
3260  PRINT  LIN(6* 

3270  SUBEXIT 
32SQ  SUBEND 

340O  SUB  In*  un- M.  Epa .  X"'*> .  Y' *  i ,  In*  >  1  Cubic  Splint  Integrator 
3480  E aed*  a* *  N< * 3/  OP  <X<1  ='  »lo>  OR  tEpi.*0> 

3490  IP  B  addt  a*0  THEN  3600 

3520  PRINT  LIN< 2 i , "ERROR  IN  SUBPROGRAM  Intun." 

3530  PRINT  USING  3540 ; N . Ep j , X *  1  * , X v H > 

3540  IMAGE  "N*  ,  DDD  ,  5a,  "Epi  *  "  ,  MZ.  6DE  ,  ",  "X*.  1  •»" ,  MZ.  6DE ,  5  X,  "X.i'N?  *  "  ,  MZ.  6DE, 

3550  PAUSE 
3560  GOTO  3480 
360O  OPTION  BASE  1 
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3610  DIM  B'.2:N-1  >,S<N>  ,G*  iSN-i  ; 

3620  FO R  1*2  TO  N-  1 
3630  H* .«  1  +  1 >-X< I-l> 

3640  ::*x<  2  >-x<  i-i ) 

3650  *•:  I/*.5*X'H 

3660  T *  <  •.  Y I  ♦  1  > - Y ^  I  >  >  •  X  •  I  ♦  1 -X  •.  I  >  *  -  \  Y  •.  I  >-Y<  1-1  >  j /'A  '■  H 

3670  S<X>«2*T 

3630  C'  I  ■*  =  3*T 

•630  NEXT  I 

3700  1 *StN>*0 

3720  w*8-4*SQR<3> 

3740  U*8 

3760  FOR  1*2  TO  H-l 

3770  T*W*<-S<I>-B<  I  3*S-:  I-l>-<.5-B<I>>4S<  I*M  >+G<  I  >> 

37S0  H*RESvT> 

3300  IF  H>U  THEN  U=H 

3620  I  >*S< I >+T 

3830  flEisT  I 

3850  IF  U  =  Eps  THEN  3740 
3860  I nt  *0 
3920  FOR  1=1  TO  H-l 
3930  H*X< I +1 >-Xi I > 

3940  Int«lm+':.5*H**YU ‘^Yvl  +  l.)  *-l-  24*H  3*<3« 

3950  NEXT  I 
3960  SUBEND 

4008  SUB  Mn*tiC2<.X2v*S'Y2'.*;,N,D> 

4010  OPTION  BRSE  1 
402O  PRINTER  IS  0 

404  0  DIM  X 4  ■ :  N  > ,  Y 4 <  M J> .  flv  4 :  N >  ,  F r  *c  4 :  N  ; ,  B  a <  5 :  N  >  t  D e rM  v  < 5 :  N  : ,  L o g d <  N  > 

4090  Ep*  i  1  on*lE-6  1  Establish  tol  Ur4nc«  rc r  convergence 
4110  FOR  1*1  TO  3 
4120  ::4<I  *.<2<I> 

4130  V4t I J»Y2< I )  !  Begin  generating  partial  area  array*. 

4140  NEXT  I 
4150  FOR  J*4  TO  N 

4160  PEDIM  X4  <  J  > ,  Y4  <  J>  ’ Incremental ly  generate  partial  area  arras'f  . 

4  170  X  4  ^  J  *  *  X  2  <  J  ' 

4  180  V  4 v  .1  >*  .’2<  J  ) 

4190  CRLL  Intun< J,Epsi lon,X4x* Y4<*>.R<J>> 

4200  Frac  D  ‘Calculate  fractional  areas. 

4210  DISP  "Generating  Rreas" , J, R< J> , Fr ac < J> 

4215  NEXT  J 

4220  FOR  K*5  TO  N 

4230  Da<k>*A<K>-ft<K>l ) 

4240  Der  i  v'IO*Da' kV<X2'-  K  >-X2'  K-l  >> 

4250  NEXT  K 

4270  FRINT  "TIME<s>  RE £R  FRRCT  R  RPER  INCREMENT  PRT£  QF  INCRERS 
E  Pf" 

4280  PRINT  LIN<2> 

4290  FOR  1*5  TO  H 

4300  PRINT  USING  4320; X2 <1 > , R<I ' , Fr ac U > , Da<I > , Der i v< I > , I 
4310  NEXT  I 

4320  I  URGE  DLDDXXXXX , DDD3XXXXX «  Z . DDDXXXXXX, DDZ. DDDXXXXXXXXXX, DDZ. DDDXXXXXNXXXXX 
.  ODD 

4330  CALL  PI otvFrac v*>,Der> wte ■ ,  M> 

4340  MRT  Y2*Y2+ < 1 E-6  >  1  Insure  nonzero  arguments 

4350  MRT  Logd*L0G<Y2> 

4365  CRLL  Pr i nt er  •: X2< *  * •  v2 •  * > . Fr ac  •"  * > » Logo <  *  > ,  N  > 

4370  CRLL  P 1  ot  2  •'  X2 ^ .  Fr  ac  V  N  ) 

4375  PRINT  PRGE 

4330  CRLL  Plot  3<X2' ♦ >,L*gd< *>,U> 

4398  PRINT  PAGE 

440O  CRLL  Plot4<Frac<4>,Ligd*'e  Ml) 

4490  SUBEXIT 
4500  SUBEND 

5000  :UB  L  t  near  <  Xt*'* ,  Y<  *  > .  M,  R.  B,  Pegs  i  ,  Re  *sa,  Tot  al  is,  Re gms, Reams , F, Bfreg.  B*V*j  t  Bf 
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tot  ,  Soon  ,  Vtr  i  tnct  ’• 


5010 

«  **■*  model:  **+ 

502O 

ON  ERROR  GOTO  Bomb 

5030 

Abort-  *0 

5040 

Y 1 *X  t  *Z*X2»Y2»0 

5050 

FOR  1*1  TO  N 

5060 

X1«X1+X<I> 

!  Sum  of 

X's 

5070 

y  i  *  v  i  ♦  r  <  I  > 

!  Sum  of 

Yf 

5030 

X2*X2+X<I>*X<n 

!  Sum  of 

X  squirts 

5090 

Y2*Y2*Yf I >4Y< I > 

•  Sum  of 

Y  squirts 

5100 

2«2*X< I >*Y< I  • 

*  Sum  of 

XY's 

5110 

NEXT  I 

5120 

Xi«Xl  N 

•  Ne*n  X 

5130 

Y1»YI'N 

!  Mean  Y 

5140 

B*<2-H*Xi*f l  >^<  X2-N*X1#X1  > 

5150 

*«Y1-B*X1 

5160 

Tot  al ss*Y2-N*Y 1 *Y 1 

•  Total  Sum  of  Square* 

5170  Rtgss«<:2-N*Xl*Yl>-2.'--X2-N»Xl*Xl.> 1  Rtgrtssian  Sum  of  Squirts 
5180  Rt*ss»Tot it ss-Rtgss  1  Rtstdusl  Sum  of  Squirts 

5190  RigmsaRtgss  1  Rtgrcssion  ftsr.  Squirts 

52O0  Rssjws*Rtiss '<n-2>  >  Rtsiduil  Mtin  Squirts 

5205  Vi r  t  inct*SORvRt»nis  -  '  X2-N*X1  *X  l  >  > 

5210  OFF  ERROR 
5220  DEFAULT  ON 

5230  F«Rtgn,s- Rtsms  •  F  Ritio 

5240  DEFAULT  OFF 

5250  Dfrtg«l  1  Dtgrtts  of  Frttdoi 

5260  Dfrts*N-2 

5270  Dftoi«Df rtg*DfrtS 

5230  SUBEXIT 

5230  Bomb:  Abort»l 

5300  SUBEND 

5310  SUB  Plotl«:X>*>,V<*),N; 

5320  1  Plot  friction  rtscttd  ritt  of  rticttor.. 

5330  ;:*»n»FNHix<X<*>,N> 

5340  X«nn»FNM:n<X<.»-,N, 

5350  Ym4**FNNix  CX<  4  > , N> 

5360  Ymn«FNMini.X(*>,N> 

5370  L INPUT  “ENTER  PLOTTER  C0tE:A«CRT,B*9872A’.Plotttr* 

5330  IF  FlotttrSU,  13«"A"  THEN  PLOTTER  IS  ’GRAPHICS" 

5390  IF  Mo«ttrfCl,n**B“  THEN  PLOTTER  IS  7, 5,  ,,98~2A" 

5400  GRAPHICS 

5410  LOCATE  5,90, 10, 75 

5420  SCALE  7mi n, Xmix, Ym i n,  Ymsx 

5430  AYES  tX» ix-Xm  n>r4,  v  Yft,i  c-Ymi  n>  4  ,  Xi»>  i  n,  Yol  n 

5440  FOR  I«1  TO  N 
5450  PLOT  X< I > , Y< I > 

5460  NEXT  I 
5470  LCPG  6 

5480  FOR  J»Xmn  TO  X«i:<  STEP  tXn>sx-X»in>r4 
5490  MOVE  J,  Ymn-Y»inr20 
5500  LABEL  USING  "D.DDE’JJ 
5510  NEXT  ) 

5520  lORC  9 

5530  FOR  t  *Vm  i  n  +  t  Y«ix-Y*1  n>  r4  TO  Y«sx  STEP  <  Ymm-Yiii  1  n  1  -'4 
5540  MOVE  Am  n-X'm  n'20,  K 
5550  LABEL  USING  “D. DDE • ; K 
5560  NEXT  t 

5570  HOVE  Xml  n .  Y«1  n-<  Yiix-Y#  1  n  >.'S 
5530  LORG  3 

5590  LABEL  "Plot  of  friction  rticttd  us  ritt" 

5595  PAUSE 

5600  IF  PlotttrBCl, 11»"A“  THEN  DUMP  GRAPHICS 
5605  IF  Plot t tr*C 1 , 1 ]«”B"  THEN  OUTPUT  7, 5; "SP0 • IN" 

5610  EXIT  GRAPHICS 
5615  GC  LEAR 
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Si  16  SUBENIT 
!6ia  SUBEND 

3630  SUB  P 1  ot  2  <  X  ■  *  '  ,  V  <■  *  '  ,  H  • 

5640  i  Plat  l-fr»C'tan  rt«t»d  1  vt  t  tms 

5630  LlMPUT  “ENTER  PLOTTeP  CODE : A-CRT , B-9S72A" , F I  at • «r * 

5660  :c  Plott«r*ri,lJ»“A"  THEN  PLOTTER  IS  “GRAPHICS" 

5670  IP  F  1  otter  PC  1 ,  THEN  PLOTTER  IS  7. 3,  “9872A" 

3630  !-  P 1 o*  t  er#I 1 ,  1 !•“ A“  THEN  PRINT  PAGE 

5690  GRAPHICS 

570O  LOCATE  10,90,10,73 

5710  SCALE  Xv5>,X<N>,0, 1 

5710  AXES  <  X(.N/-:<v5j  W4,  .S3, OS'1, 0 

3738  FOR  I«3  TO  N 

5740  PLOT  Xt  I  )  ,  1 -Y f l j 

5750  NE:  T  I 

3760  lORC  6 

3770  FOP  J«X'-3>  TO  X*. N  •  STEP  ■'N  N)-X<5<  >r4 

3780  HOVE  J.-.OI 

3730  wABEL  USING  "DODD" i  J 

5300  NEXT  J 

5010  MOVE  (XtN'-XvS)/  2*X'  3  ■, 07 
33i0  LABEL  "Tibi  <*eec.nd*3” 

3830  lORC  8 

5840  FOP  K.,25  TO  1  STEP  .23 
3850  MOVE  •  X'  N>-:i <5>  >'  20*X''5> ,  K 
3860  lABEL  USING  “2.DD“;K 
3870  NEXT  K 
5838  LORC  3 
3890  lDIR  F I '2 

5900  MOVE  fUNT-MS  •  'v20*X(3».  .03 
3910  LABEL  “Friction  Unrcicted" 

3913  PAUSE 

3920  IF  Plott«r*Cl,:3«"A'  THEN  I  VIP  GRAPHICS 
3923  IF  Plotter*! 1, I }>"B"  THEN  OUTPUT  7,S;“SP0rIN" 

3930  EXIT  GRAPHICS 

3931  GCLEAP 
3935  SUBEXIT 
3940  SUBEND 

5930  SOB  Plot3<X<»*,Y<*>,ti> 

5960  1  PI  :>t  o  f  log  di  ipl  tcetient  vs  time 

5970  DIM  A'.N'.BCN)  'P»rti»l  Rrrty*  for  Kinetic*  C  »1  eul  »t  i  or, 
3980  lINPUT  “ENTER  PLOTTER  COPE:  A«CRT,  B*987iA" , P I jt tf r* 
3990  IF  P1ott«r*Cl,n*"A“  THEN  PLOTTER  IS  "GRAPHICS" 

6000  IF  F I ot  t  er*t 1 ,  1  J«*S"  THEN  PLOTTEP  IS  7(5("gj72f,“ 

6010  GRAPHICS 
6020  LOCATE  10.90,10.73 
6030  SCALE  Xt5>,XtN"',0,5 
6035  lPIP  0 

6040  AXES  <X<N»-Xt5>>'4, 1,X<3),0 
6030  FOP  I«3  TO  N 
6060  DRAW  X< I > , Yt I > 

6070  NEXT  I 

6088  LOPG  6  '  Libel  ’  N«  *xe* 

6090  FOR  I-XtSi  TO  X<M)  STEF  i V N)-X<5) )r4 

6100  MOVE  I ,-.03 

6110  LABEL  USING  "DDD“;I 

6120  NEXT  I 

6130  HOVE  (X(M>-X<3>>--2»X<3',-.3 
6140  LABEL  "Tim*  <*«c>" 

614*  LORG  8 

6150  FOR  7* 1  TO  5 

6160  MOVE  Xi5>-.01*fX<M>-X'5‘>,J 

6170  LABEL  USING  "D.D“;J 

6  ISO  NE.-.T  J 

6139  LDIP  PI.  2 

62O0  LOPG  6 
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6216  MOVE  X<5>*.05#<X<N>-X<3>>,2.3 
€229  LABEL  'Lntdi »pl 4c«»*n«  >" 

6230  LDIR  9 
6240  PAUSE 

6250  IF  P1ott«r4Cl, 13»"A“  THEN  DUMP  GRAPHICS 
6255  IF  Plott«r#tl, THEN  OUTPUT  7, 5; “SPO- IN” 

6260  EXI T  GRAPHICS 
6265  GCLEAR 
6270  PRINTER  IS  0 

6275  INPUT  “ENTER  TWO  POINT  NUMBERS  FOR  THE  SLOPE  CALCULATION" , S, T 

6280  Nptid-T-S»l 

6285  REDIM  Av 1 :Nptsd», Bt ltNptid) 

6280  PRINT  LIN-33 
6285  FOR  1-1  TO  hptid 
6300  A(  I ,'«X< S*I~1  > 

6305  B< I>»V<S*I-1)  . 

6310  NEXT  I 

6320  CALL  L>n«4r < At*) , 8<  *  '  ,Npt id, Int , Stp,C0,Ct ,  C2,  C  3,  C4 ,  C5,  C6,  C7,  C8 ,  C9,  Var  > 

6330  PRINT  USING  6340;  A <  1)  ,  AOlpt »d  ) .  S  1  p,  Var 

6340  IMAGE  “Slop*  of  plot  between  times  ".DDDD,"  and  “,DDDD,"  sa.ondi  it  "-'D-DDE 
".,  D.  DDE ,  “  14C-1" 

6350  PRINT  LIN<3> 

6360  LINPUT  “DO  YOU  WANT  TO  REPEAT  THE  CALCULATION  USING  OTHER  POINTS?" , Rp* t 
6370  IF  (Rptftl, 13»“Y">  OR  <Rpt#tl, 13«“u“)  THEN  GOTO  6275 
6375  SUBEXIT 
6360  SUBEND 

6600  SUB  Plot<X<4>, Y<*>,N> 

6610  REM  Plot  dAxdt  v*  A  and  chooia  linaar  portion  for  kinatic* 

6623  PRINT  PAGE 
6625  OPTION  BASE  1 
6630  PRINT  LIN<3> 

6640  PRINT  "Kinttici  plot  of  dx'dt  w»  x.  x  »  friction  daconpoaad" 

6650  DIM  Xl<N',VHN) 

6660  C*FHM*v.<Y<«>,N,5) 

6680  F*C  10 

6700  PLOTTER  IS  "GRAPHICS" 

6710  GRAPHICS 
6720  LOCATE  10,90,10,75 
6730  SCALE  0, I.O.C 
6740  AXES  .25,2.5*F,0,0 
6750  FOR  I«5  TO  N 
6760  MOVE  X< I ) , Y< I J 
6770  LABEL 
6780  NEXT  I 
6780  LOPG  6 

6800  FOR  1-0  TO  1  STEP  .25 
6910  MOVE  1,-.»*F 
6320  LABEL  USING  "Z.DD";l 
6830  NEXT  I 
6840  LORG  8 

6850  FOR  1-2. 5«F  TO  C  STEP  2.5*F 
6360  MOVE  -.01,1 
6870  LABEL  USING  "02.09":! 

6888  NEXT  I 
6380  LORG  6 
6800  MOVE  . 5, - . 8»F 
6810  LABEL  "Fraction  Dtcgi*po»«d" 

6920  LDIR  PI -2 
6930  MOVE  .025,D+3»F 
6940  LABEL  "Palativa  Rata" 

6950  LDIR  0 
6960  PAUSE 
6970  DUMP  GRAPHICS 
6900  CCLEAR 
6990  EXIT  GRAPHICS 
7000  PRINT  PACE 
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7310  PRINT  LIIK2) 

7010  FPINT  "Lin«»r  R«gr«i*icn  Ar»lv*i«  of  th*  Dati" 

7033  PRINT  LIN<2; 

7040  INPUT  “ENTER  BEGINNING  ANT  ENDING  POINT  FO®  CURVE  FITTING'',  B*gt  n, End 

7050  - t **End-8tg t n* 1 

7060  PEDIM  XUlsPUi.VliliP'il 

’070  -OP  1*1  TO  Pts 

7000  I'  I  Eegin»I-l.' 

7090  VI  >.  I  ■» V<  5*gi  n-*I-l  > 

7100  NEXT  I 

7110  PRINT  USING  7120;Xl.l),XHPt*> 

7120  IMAGE  “Plot  ij  b«tw<*n  tiMi  ".DDDD.D,"  s#t ondi  and  '', DDDD.D,''  ascond*. 
7130  PRINT  L I N < 2 > 

7140  CALL  LiniirtXl  *) , Y1 <*> , Pt *, A. B, C0, Cl , C2, C3, C4, C5. C6,C7, C8, C9, V> 

7150  PRINT  USING  TlSOjABS'EV.V 

7160  IMAGE  “Rat«  Comtant  it  “,D.DDE,  “**c-l.  M«an  Variance  1*  ”,D.DDE 
7170  PRINT  1_  I N  <  3  > 

7130  X»ix*FNHixOa>.*',Pti,  1  > 

7190  <  VI  (  «  >  ,Pt»,  l  ’ 

7200  Xai  n*FNMl  n.  XI  v  *  1 ,  P*  *> 

7210  V«n  n*FNMi  ru  VI  <*  > ,  Pt  a  • 

7220  E* > Xii-i-XiM n>  •  10 
7230  F*>Y»ix-Yii.in'  10 
7240  PLOTTER  IS  “GRAPHICS" 

7250  GRAPHICS 

7260  LOCATE  20,100,20,30 

7270  SCALE  X»i n, Xnax, Y»i n,  Yaay 

7200  AXES  2*E,2*F,X«.in,  V«in 

7290  FOR  1*1  TO  Ft  4 

730O  MOVE  XUD.Ylfli 

7310  LABEL  “♦* 

7320  NEXT  I 

7330  FOR  J»1  TO  Pt* 

7340  Yl<  J>»X10>*BfA 

7350  NEXT  J 

7360  MOVE  XI < n ,Y1 '  i> 

7373  FOR  1*1  TO  Pt i 
7330  DRAW  XU  I  >,V1'  I 
7390  NEXT  I 
’400  ,-OPG  6 

74  10  FOR  I*Xmir,  TO  Xu, ax  STEP  2*E 
7420  MOVE  I,Ymin-.l*F 
7430  LABEL  USING  *2. ODD"; I 
“440  NEXT  I 
7450  LORG  8 

7460  FOR  !*Yl»m*2*F  tq  ,«,»*  STEP  2-F 
7470  MOVE  X« 1 n-. G4-E , I 
7480  LABEL  USING  "DZ.DD":I 
7490  NEXT  I 
7500  LORG  6 

7510  HOVE  X*in*5*E,v»in-P 
7520  LABEL  "Friction  D«coi»poa*d" 

7530  LDIR  PI-2 

7540  MOVE  X*nr.».4»E,Y»in*5*F 

7550  LABEL  ”R«latu>«  Rata" 

7560  »AUSE 
7570  DUMP  GRAPHICS 
7580  GCLEAR 
7590  EXIT  GRAPHICS 

7600  LINPUT  “DC  YOU  WANT  TO  RErE AT  THE  CALCULATION’", Pop# it* 

7610  IF  R*p«it<C  1 .  1  THEN  GOTO  7000 

7620  SUSEXIT 

’630  SUBEND 

7700  BEF  FMMonf Av*> , N> 

77ie  Ml n*A< 1 > 

7720  FOR  1*1  TO  N 
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7730  IF  PtIXHin  THEN  Hin-P<I> 

7748  NEXT  I 
7750  RETURN  Min 
77E0  FNEND 

7770  SUE  Pr  i  nt*r  ■ Ti net * ) ,  Zi <  spt  *> ,  Frac  v* ' ,  LogJC*^ ,  N> 

77S0  !  Print  out  viiuii  used  in  kinetics  cklcuUtiont 

7790  DI3P  “CHOOSE  PRINTER  PND  PRESS  CONT” 

7800  PPUSE 

7310  PRINT  "Pvt  Tinieitec/  V  LoaO'l  a  <1-sj  L 

n  lQO«(l-a'“  ~  - - 

7420  PRINT 

7830  FOR  I«5  TO  N 

7835  IF  Fr»c<l>-1  THEN  GOTO  7870 

7840  PRINT  USING  7860;  I ,  T  i  me  <  IT ,  Di  ip  <  I ) ,  Logd(  1 > ,  Free  <  I ) ,  1-Frac  <  I) .  LOG<  100*<  1  -Ft  a 

c<D') 

7850  NEXT  I 

7S60  IMPGE  DDDXXXXX, BDDDXXXXXX, DIB. BXXXXX.MZ. DDDXXXXX, Z. BDDXXX,Z. DDOXXXXX.MZ.  D 

DD 

7370  PRINT  PPGE 
7830  SUBEXIT 
7890  SUEEND 

7900  SUB  Plot4<X<*>,Y<e),Z.’ 

7901  DIM  U<Z),V(2)!  order  plot  of  1 og< di sp I acement )  vs  1 ogt 100*< 1 -a> > 

7905  LINPUT  “ENTER  PLOTTER  CODE:  CRT»R,  9872P-E" , PI ot t #r* 

7910  IF  PI  otterfCl ,  U»“fi“  THEN  PLOTTER  IS  “CRPPHICS" 

7915  IF  Plotter SCI, 1 J-“B"  THEN  FLOTTER  IS  7,5,“9872P" 

7920  MPT  1  )-x 

7925  MPT  X-<10S)#X  !  Eliminate  negative  logarithms 

7930  FOR  I *5  TO  Z 

7933  IF  X<I>«0  THEN  GOTO  7937 

7935  X< 1 t-L0G<  Xt I > ) 

7937  NEXT  I 
7940  P-5 

7950  C-5 

7970  LOCPTE  15,90,10,80 
7980  SCPLE  0.0, 0,C 
7905  SHOW  0,p,0,C 
7990  PXES  1. 1,0,0 
8080  FOR  1-5  TO  Z 

8092  IF  Y  <  I  )<  0  THEN  Y<I)-0  i  Ignore  negative  1  ogar  i  t  hits . 

8803  IF  :<<I)<0  THEN  XCIT-0 
8010  PLOT  XU>,Y<t> 

8020  NEXT  I 

8825  LORG  6 

3030  FOR  1*0  TO  fl 

8040  MOVE  I,-.02*C 

8050  LPBEL  USING  "Z.D";I 

8060  NEXT  I 

3070  MOVE  P'2, 08<C 

8080  LPBEL  “Lnt 108» t 1 -a) ) " 

8085  LORG  8 

8098  FOR  I-l  TO  5 

8108  MOVE  - . 32*P, 1 

8110  LPBEL  USING  "D.DD";I 

8120  NEXT  I 

8130  LOIR  PI  2 

8148  MOVE  . 0S*P, . 5*C 

8150  LORG  6 

8160  LPBEL  "Lntdi sp 1 ac stent > " 

8170  PPUSE 

8180  IF  Plotter*! 1, 1 3-“P"  THEN  DUMP  CRPPHICS 
8190  IF  Plotter*Cl,13-"8“  THEN  OUTPUT  7, 5;  ”  SPO-'  IN" 

8200  LDIR  0 
8210  EXIT  CRPPHICS 
3220  GCLF*P 

3230  INPIP  “ENTER  TMO  POINTS  FOR  SLOPE  CPLCULPT ION" , B, D 
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3343 

3353  »El!n  u ' i : Opt i  ' , v ■ i : op*  s  > 

8060  rOK  I-l  TO  Op<  » 

8070  0-  I  ■«Xt»*I-I) 

8088  V  1  ■  »y <  B* I - 1 > 

8288  NE"T  I 

8708  GALL  Llt.t  v  U'  Cot*.  Int ,  SI  op«, C8, C  1 ,  C2 ,  C  0.  C •*.  C5,  C 

3318  PS'NT  L1M<3> 

8308  CRINT  USING  8333;  U‘  1  '• ,  Ut  Opt  4>  ,  SI  opt.  V*r 

8338  IMAGE  "Slop#  of  order  curve  petueen  vtljti  "^D.DD."  end 
■  ♦  -  ■.r.ni'D 

8343  .I>.?UT  '30  YOU  NANI  T0  PEPEAT  THE  C  AuCUlA  t  :  **  .  A  ;  e  •  r  f 

3350  I r  Agei  nJI  1 ,  1  ]>  *Y"  THEN  GOTO  8*033 
3363  -UBS  IT 
8378  SUBEND 


.C7,C3 
,9. DD, 
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APPENDIX  C 

Program  for  Acquisition  and  Analysis  of  Thei*ogravi*etric  Data 
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FPOCPA'i  FOF  AQUISITION  AND  ANALYSIS  OF  CONSTANT  HEATING  RATE 
DIFFERENTIAL  SCANNING  CALORIMETER  DATA. 


10  OF ’IO*i  EASE  1 

10  S!i  ■  •  lO00',Te».„'  1080  ,:.i<  1000:  ,V3<  1000  > 

30  FRINTEP  IS  IS 
40  PR 1 h T  LlNc4> 

50  PRINT  "Data  aquisitior.  program  for  tht  Sc  »nr  i  Di  f  ftrtnt  I  *1  * 

SO  FSINT  "Caloriattar.  Ust  ictnmnf  digital  voltmeter  lint" 

70  FAINT  “9.  Ordinate  a  is  t*  calibrated  for  10  millivolts" 

30  PAINT  “full  seal*.  Tiaptrtturi  is  calculi'  tl  si  TIME  “ 

90  FAINT  'HEATING  FATE.  Data  is  taken  *•  *  hi  appro  i  ...  st  4  ntr1 
100  oc  INT  “of  or#  po'nt  ptr  4  i,con#i.  Edit  line  310  ‘o  change  rate" 

110  PFIMT 

123  FAINT  “Follow  t  ht  CAT  prompt*.  Prtss  CONT  *  c  #r*e>  data* 

130  FAINT  ‘into  th#  computer  or  to  have  s  plot  or  "he  OFT  put" 

140  PAINT  "onto  t  h#  thermal  pnntir.  Ttniniti  t  F.«  pregr  is  b"  pressing" 

145  FAInT  "STOP.  Edit  K5  if  you  want  additional  graphics  f«i*ur#i,N 

146  FAINT 

150  PAINT  “EDIT  THE  FAOCRAH  AND  PRESS  CONT" 

155  PAUSE 

ISO  FAINT  PAGE 

170  FAINTER  IS  0 

ISO  DIN  AFC603, EJIE03.CSC8OJ 

190  LINPUT  "  ENTER  TEST  IDENTIFIER,  LINE  l  “ ,  A* 

200  L INPUT  “ENTER  TEST  IDENTIFIER.  LINE  2“.E* 

210  INPUT  ENTER  INITIAL  TENFEFATUAE  OF  THE  RUN  v  ".Tinn 
210  INPUT  "ENTER  THE  HEATING  FATE  IN  DEGFEES  C  ptr  MINUTE", R*t# 

225  INPUT  “ENTER  NASS  OF  SAMPLE  cst'*,SMpUs*(i 

226  INPUT  "ENTER  SENSITIVITY  >  He  *1 -sees  * ,  S«ns  i  1 1  vi  t  y 

230  A»t«»Rs»«  €0  'Convtrt  rst#  »o  d#gr##s  p»r  second. 

240  PRINT  L I M <.  3  * 

250  °RINT  At 
260  PRINT  Bt 

270  "PINT  USING  230;T,nit,Rv«»£0 

275  PRINT  USING  235; Sapp I email , Sensi t iui t y 

2S0  INAGE  'Imtisl  T empe "iturf  is  “,DDD,"  K  Heati'g  R»t#  is  “,DD,"  degmiin" 
285  INAGE  "lias*  of  JMipIt  is  ".DD.DDD,"  mg.  Sensitivity  is  “,DD,"  meal  s»c " 
290  PRINT  L INC 2 * 

300  REN  Th#  following  lints  cor.figurt  tht  sc  inning  digit*)  voltmeter  *nd  t  i  at 
clod  and  stt  key  irttrupts. 

310  OUTPUT  701; “R  M0  0”!Configurt  DVN  for  inttrupt  transfer. 

320  ON  ERPOR  GOSUB  Esc  apt  I En*b 1 1  run  trrors  to  bt  corrected. 

330  ON  PEY  *7  GOSUB  End'Operator  ttraination  of  data  aquitition. 

335  ON  KEY  «8  GOSUB  Exit 

340  ON  INT  #7  GOSUE  Read 1 Branchi ng  routine  for  data  transfer. 

330  OUTPUT  9; "A- U1*I2" 'Configure  clock. 

360  CONTROL  MASK  7; 1 23 'Look  for  inttrupt  flag. 

370  PRINT  L I Nt 2 ) 

373  PRINT  "TEMPERATURE  <K/  ORDINATE  PCINTR" 

380  PRINT 
383  I«1 

398  DISP  "Prtss  CONT  to  begin  data  aquisition" 

400  PAUSE 

410  OUTPUT  9; "Ulu“ I  Start  th#  clock 

420  DISP  “PRESS  k  7  TO  TERMINATE  DATA  GATHEAINC” 

430  CARD  ENABLE  7'Allow  DVN  inttrupt*  to  t*kt  pi  act. 

440  IF  !<I000  THEN  GOTO  440 
430  DISABLE 
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460  PRINTER  IS  16 

470  PRINT  "ALLOCATED  COPE  STORAGE  IS  FULL" 

480  PRINTER  IS  0 

490  GOSUB  End 

300  Eiit:  OUTPUT  701;  •«" 

303  STOP 

310  Raad:  1  INPUT  DATA  FROM  BSC 

310  ENTER  701 JY' I' 

323  Y  .  I  )»Y1' l  >*1000*20  ‘Convart  milliuoltt  to  »  0  to  !0  tcala. 

330  OUTPUT  9; "U1V" 

540  ENTER  9; Ti ■* 

330  Tamp>  I^*Ti««*P»t*1000*TinH 

560  PRINT  USING  370; Tamp<I > , Y;  l > , I 

370  image  ddddxxxxxxxxxxxxxxx,mbd.ddxxxxxxxx,dddd 

580  I • I » 1 

630  OUTPUT  701 ; “0“  . 

660  CARD  ENABLE  7 
670  RETURN 

680  End;  !  Procatt  Bata 

683  '■' <  1  > «Y <2 >  (Bad  dill  point  aiisint'or 

690  Tajnpi  l  )«Ti  n  1 1  'Mata  iur<  firtt  point  it  corrtet 

708  Nptt*!-2  'Etiminata  axtranaout  data  point! 

710  REBIM  YCNpt t > , T»mp i Npt t >  'Cond«nt«  tha  trrtut. 

720  Mi  ni  mu»«Fnm  i  n  c  YC  *  > ,  Npt  t  >  'Calculata  louet'  '.s’taga  output  from  DSC. 

730  MAT  Y«Y-tMi  o'  mutt)  'S*t  zar  :  point  tt  louut  output  uoltaga 

770  CALL  CrapA'A*.  B«,  Taap<*>,  Yt  *'>  ,Npt  t> 

790  LINPUT  “DO  YOU  MANT  THE  DATA  STORED?’, C* 

790  IF  C#Cl,ll»"Y"  THEN  CALL  St  or* ( A*, B* , Ta*p< »  , Y< *> , Npt * , T i ni t , Rat ♦, Sampl m 
tt ,  S«nti  t  i  ui  tv*> 

80O  LINPUT  "DO  YOU  WANT  TO  EXPAND  THE  FLOT",D» 

310  IF  DBtl.lJ-’Y"  THEN  CALL  Expand'’  A*.  B*.  Taap.  *) ,  Y<»  >  > 

828  LINPUT  "DO  YOU  WANT  TO  MAKE  A  LINEAR  BASELINE  CORRECTION"*"  ,E* 

330  IF  EIC1,  n«"Y"  THEN  CALL  Tr ant f or mC A# . B», Tampt* > , Yi *> , Npt t > 

833  PRINTER  IS  16 

836  PRINT  PAGE 

837  PRINT  "Paak  intagrator  uting  »n  arbitrary  bataltna  correction" 

838  PRINT 

839  PRINT  "3***  point  nu«b«rt  corrttpond  to  origional  d»t »  lilt" 

840  PRINT  LIN<2> 

841  PRINT  "TM*  btitl  mt  i  i  calculated  using  a  ttaight  line  f  rot  the  firtt  to 
i»i»  1  »st  data  point" 

842  PRINT  L  l  N  <  3  > 

843  PRINT  "Pratt  STOP  if  vou  witK.  to  terminate  the  program.  Pratt  CONT  if 
you  with  to  continue." 

844  PAUSE 

845  PPINT  PAGE 

848  INPUT  "ENTER  BEGINNING  AND  ENDING  POINTS  FOR  NUMERICAL  INTEGRATION" , Bagi n, 
End 

038  Opt t*End-B* gi n*l 

032  PRINT  USING  034;Qplf 

834  IMAGE  "ENTER  1  AND  * , ODD, "  FOR  INTEGRATION  LIMITS* 

836  REDIH  X3C Opt t > , Y3< Opt t > 

838  FOR  L* I  TO  Qptt 

860  Y3'L>«YiL*B«gin-l > 

862  x3iL>*Tamp<L*Bagin-l > 

864  NEXT  L 

866  S 1  ope3«< Y3<0pt  a  >-Y3t  1  •  >  •  t  >:3  ( Opt  t  >-X3  <1 )  > 

867  B*Y3<1> 

866  FOR  L«1  TO  Qptt 

870  Y3<L>mY3<L.'-Slopa3»<X3<L'-'0<  I>»-B 

072  NEXT  L 

874  C*«B*S":  Expandad  plot  peak  intagration" 

876  CALL  Graph2*  flf,  C*,  XJ<  *>  ,y3<.»>  , Qptt,  Area) 

875  LINPUT  "DO  YOU  WANT  TO  CONTINUE  WITH  K INET ICS?" . St  op* 

830  IF  Slop*C 1 ,  li<  >"Y*  THEM  COSUB  Exit 

881  CALL  M  net  i  cf  R»,C»,  X3>  *  ,  V3r»  > ,  Op 1 1,  Area' 
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304  GOTO  300 

363  EicapPl  1  iorric;  arv  run  irrori, 

660  PRINTER  IS  16 

300  PRINT  t«n«r  0EiEE0TEI.  EDIT  PROGRAM  AND  PRESS  COST." 

310  “ROSE 

920  PR INTER  IS  0 

930  RETURN 

1002  SUE  St  or«  ■  A» .  tl  N.  TO,  R,  S«>,  S«r,s  ' 

1010  OPTION  IPSE  1  'Suprjuttns  for  storing  da*  a  (r;n  X-Y  plotter  input. 

1013  SHORT  R • N  • ,  S  \  N  » 

1020  lINPUT  “ENTER  MASS  STORAGE  MEDIUM" , F# 

1630  ..input  -enter  6  character  file  name" , g* 

1040  G*»G*S.“:“«.Fi 

1030  CREATE  C*,Mr9 
106O  ASSIGN  42  TO  G* 

1070  PRINT  42;  fit,  SI.  N,  T0,  R,  SN,S»ns 

1075  MAT  R«X 

1076  NAT  S»Y 

1080  MAT  PRINT  RUF.S 

109O  PRINT  LIN<3> 

1100  PRINT  "DATA  IS  STORED  IN  FILE  “LCHR*< 34 ' ; G»; CHR4 > 34) 

1 1 10  ASSIGN  42  TO  ♦ 

1120  SUBEXIT 

1130  SU8END 

1200  SUB  Graph' Af,B4,X'*),Y<*',N> 

1201  OPTION  BASE  1 

1203  INPUT  "ENTER  PLOTTER  CODE:  1-CRT;  2-9672A” , Coda 
1206  IF  Coda«2  THEN  GOTO  1270 

1210  PRINT  PAGE  'Su&rcutina  for  plotting  data  from  an  X-Y  plottar 

1220  PRINT  LINO  ’  'using  th*  mtarnal  tharaal  prlntar. 

1240  PPINT  A* 

1230  PPINT  B* 

1260  PRINT  l:n‘3> 

1270  DIM  CSI661.DSC601 

127!  LINPUT  "ENTER  X  AXIS  LABEL  IF  OTHER  THAN  TEMPERATURE  <IO",Ct 

1276  IF  C4«““  THEN  CS*" TEMPERATURE  <«)* 

1277  LINPUT  "ENTER  Y  AXIS  LABEL  IF  DESIRED", D4 
1276  IF  t»-"’  THEN  D*-“Ra1ativa  Millivolts" 

1280  F*FNMat<Y<«i,N) 

1290  B*XfN>-X<n 
1300  A*F 

1310  LINPUT  “DO  YOU  WANT  A  MRXIMUM  Y  USED  FOR  THE  GRAPH  INSTEAD  OF  THE  CAlCULAT 
ED  MAXIMUM", J4 

1320  IF  J4H,1Jp*Y"  THEN  INPUT  “ENTER  i«sx',H 

1339  IP  J4C  1 , I 3«"Y"  THEN  A»H 

1340  IF  Cod««2  THEN  PLOTTER  IS  7,3,"9872A* 

1330  IF  Cod4-l  THEN  PLOTTER  IS  "GRAPHICS" 

1360  GRAPHICS 

1370  LOCATE  15,90,13,73 
1380  SCALE  X<1>,X<N),0,A 
1390  AXES  Br4,Arie,X<l>,0 
1400  HOVE  X'1>,Y<1> 

1410  FOR  I»1  TO  N 
1420  DRAI!  X<  I  >  ,  Y<  I  > 

1438  NEXT  I 
1440  PENUP 
1480  LORG  $ 

1490  FOR  1«X<1>  TO  X<N>  STEP  i-  4 

1380  MOVE  I , -Ar 100 

1318  LABEL  USING  1330; I 

1320  NEXT  l 

1330  IMAGE  ODD 

1340  LORG  0 

1330  FOR  I»Arl0  TO  A  STEP  A'10 
1360  MOVE  -»rl00*X' 1),I 
1370  LABEL  USING  1590; I 
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1560  NEXT  ’ 

1590  IMAffi  DD.  D 
1000  LORG  6 

1010  MOVE  B'2*X<;i>,-A  15 
1020  LABEL  TRIHFcC*' 

1630  LDIR  P I '2 
1640  L  0  PC  4 

1650  MOVE  2 

1660  LABEL  TPIM*<D*> 

1670  PAUSE 

1680  IF  Cod**  1  THEN  DUMP  GRAPHICS 

1690  EXIT  GRAPHICS 

1700  GCLEAR 

1710  SUBEXIT 

1720  SUBEND 

1730  DEF  FNMax<A(*»,N) 

1740  Max*A< 1 ) 

1750  FOP  1*2  TO  N 

1760  IF  A <  I  )  >M*.x  THEN  Max*A<I> 

1770  NEXT  I 
17S0  PE.TURN  Max 
1790  FNEND 

2OO0  SUB  Expand<C*.D*,X<*j,Y<:*>> 

2010  OPTION  BASE  1 

2020  DIM  Z<1O00>,T<’ ieO8>,RdO08>,Hd000< 

2025  DIM  1*1803, RBC80I 

2030  INPUT  “ENTER  LOWEST  AND  HIGHEST  DATA  POINTS  FOP  PLOT  EXPANSION  rP 1 , P2 ' " , P 1 
,P2 

2040  N*ABSCP2-Pn*l 
2050  REDIM  2<N1,T<N) 

2060  FOP  1*1  TO  N 
2070  ;<i>*x<i*pi-i> 

2030  T< I >*Y< I *P1 -1 > 

2093  NEXT  I 

2100  I»*D*l“-EXPANDED  PLOT" 

2110  LINPUT  "IS  THIS  A  BASELINE  COPPECTED  PLOT?",N* 

2120  IF  MBCl.lJd-V"  THEN  GOTO  2175 

2138  IF  M*C1.1J*"Y*  THEN  CALL  Pr i ntz<C», l#, 2<*> , T<*> , N ■ 

2140  IF  M*Cl,tJ«"Y"  THEN  CALL  Cr  aph2<  C* ,  I  #,  Z<  *  >  ,  T<  * .) ,  M,  Area) 

2141  PRINTER  IS  16 

2143  PRINT  “DO  YOU  WANT  TO  INTEGRATE  THIS  EXPANDED  PLOT  USING  ANOTHER  BASELINE 
CORRECTION''" 

2144  LINPUT  N* 

2145  PRINT  PACE 

2146  IF  Nltl.DO’Y"  THEN  GOTO  2170 

2147  INPUT  “ENTER  THE  DATA  POINT  NUMBERS  <n«w  1  i  *t  >  FOR  ORDINATE  POINTS  “,P3,P4 
2149  Mx*P4-P3*l 

2149  REDIM  R(Mx) , W<M\> 

2158  PRINT  USING  2152;Mx 

2152  IMAGE  "USE  1  AND  * , DDD, "  FOR  INTEGRATION  LIMITS" 

2153  FOR  K*1  TO  Mx 

2154  l»<K)*Z<K*P3-l> 

2155  Mil  <*TiK*P3-l> 

2156  NEXT  K 

2157  R**"Expand«d  plot:  Arbitrar’.'  Budint  correction  for  curve  integration." 
2153  S 1  op*  2*<U<Hx)-H<l>)x<R<H*;-R<l>> 

2159  lb*W< 1 > 

2162  FOP  K ■ 1  TO  Mx 

2163  H<KJ*H<K)-Slop*2*<RrK>-RU)>-Bb 

2164  NEXT  K 

2165  CALL  GrapH2<C$,R*,Rr»),M<*>,Mx. Ar*») 

2170  IF  M*C1,1I-"Y"  THEN  SUBEXIT 

2175  CALL  Graph-CB, IB, Z<*» , T<*> . N> 

2180  SUSEXIT 
2195  SUBEND 

2190  SUB  Tran*fcr*.F#,G»,X<#>,Y<*<,N> 
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2200  OPTION  BASE  1  1  Linear  baseline  correction  for  2  .‘ariable  plots. 

2210  BIN  YntwiN) 

2212  INPUT  "ENTER  POINT  NUMBER  TO  BEGIN  BASEuINE  CORRECT  ION" , Opt 
2215  MH  HRC80] 

2220  S»tYi  Opt  t+Y<Qpt +1 i  (Opt  t-2  >  V  3  ’  Average  three  points 
2230  V»<V<M >eYt M-l >+Y<N-2> >x3  i  for  s 1  ope  calculation. 

2240  V.  Ope  =  >  V-S  '•  ••  <  X  >.  N  >  -X  cl  > ) 

2250  CCP  1=1  TO  N 

2260  YneuC  1  •  =  Y  <  I  Y-S 1  ope  »'  X;  I  >-X<  1  >  ;-S 
2220  NEXT  I 


2280  HS*GS!J. "  :  L  i  near  baseline  correction" 

2290  CALL  Cr  aph2  <F»,  HI,  X<'*>  ,  Ynew<  *>  ,  N,  Area; 

2330  LINPUT  "BO  YOU  WANT  THE  PLOT  EXPANDED' N* 

2335  IF  M»tl,13="Y"  THEN  CALL  Exp  andtF* ,  H*,  X<  *>  ,  Ynewt  *>  > 

2340  SUBEXIT 
2350  SUBEND 

2400  DEF  FNMaxabs »'  A<  *>  ,  N  > 

2410  Ha.'  =ABS 1  A 1  1  >  > 

2420  FOP  1=2  TO  N 

24  jo  IF  AB3;  At  I  '  >  .'Max  THEN  M ax*ABS < A t I  . 

2440  NEXT  I 

2441  RETURN  Ma- 

2442  FNEND 

2450  DEC  FNI-Imt  Bt»  >,N  ‘ 

24C0  Ml  n  =  B  ‘  1  ■> 

2470  FOR  1=2  TO  N 

2480  IF  BvIXHin  THEN  Mtr,«B<.l> 

2490  NEXT  I 
250O  RETURN  Min 
2510  FNEND 

2600  SUB  Gr apH2i A», B4, M< ♦  ; . Y<» ' , N,  Int < gral > 

2881  OPTION  BASE  1 
2805  PRINTEP  IS  0 

2810  INPUT  'ENTER  PLOTTER  CODE:  1«CRT;  2=98720" , Code 
2815  IF  Code *2  TnEN  GOTO  2670 

2820  PRINT  PAGE  '  Subroui  int  for  plotting  data  frcin  an  N~Y  plotter 

2830  PRINT  LIN*-  3  )  'using  the  internal  thermal  printer. 

28 40  PRINT  A* 

2650  PRINT  B* 

2660  PRINT  LIN<  3; 


26-0  DIM  C IC80] , D4C801 .UxC 10001, Vyt 1000; 

2671  PEBIN  U>  < N  J , Vy  t  N > 

2675  LINPUT  "ENTER  X  AXIS  LABEL  IF  OTHER  THAN  TEMPERATURE  <K>", C* 

2676  IF  CB  =  >"  THEN  C *• " TEMPERATUPE  CIO" 

2677  LINPUT  "ENTER  v  AXIS  LABEL  IF  DESIRED". D» 

2678  IF  D*«“”  THEN  D*«"Relaliue  Millivolt*" 

2680  F =FNMa>  abs  <  Y*  *> ,  N > 

2690  B  =  X f  N  > -X  c  1 ) 

2700  A»F 

2710  LINPUT  "DO  YOU  WANT  A  MAXIMUM  Y  USED  FOR  THE  GRAPH  INSTEAD  OF  THE  CALCULAT 
ED  MAXIMUM", J* 

2720  IF  t*Cl,tl’"V  THEN  INPUT  "ENTER  Y»a*",H 

2740  IF  JfCt,lJ=*Y"  THEN  A*H 

2750  IF  Code* 1  THEN  PLOTTER  IS  "GRAPHICS" 

2755  IP  Code *2  THEN  PLOTTER  IS  7,5,n9872A" 

2760  GRAPHICS 


2770  LOCATE  15,90,15,75 
27J0  scale  Y«l>,XcM>,-fl.A 
2790  AXES  B'4,A  ie,xa>,0 
2000  MOVE  Xtl),Y<l> 

2010  FOR  1  =  1  TO  II 
2020  DRAW  XcD.YCI) 

2030  NEXT  I 
2040  FEMUR 
2880  LORG  6 

2890  FOP  ;«:<«  I;  TO  X<N>  STEP  B  4 


62 


NWC  TP  6337 


i 

I 


1*08  MOVE  1,-0- 100 
2910  LABEL  USING  2939:1 
2920  NEXT  I 

2939  IMAGE  ODD 

2940  L0RG  S 

2939  FOR  I  =-A  TO  A  STEF  A'10 
2960  MOVE  -J  100*X(I),I 
2970  LABEL  USING  2990:1 
2930  NEXT  I 
2990  IMAGE  MDD.D 

3000  LORG  6 

301O  MOVE  B.'2+XM),-A'15 
3O20  LABEL  TRIM*tC*> 

3030  LDIR  PI  2 
3040  LORG  4 

3030  MOVE  -8'10+X<l;tA  2 
3060  LABEL  TRIM*<B*> 

3962  LINPUT  "BO  YOU  WANT  TO  INTEGRATE  PART  OF  T HE  CJR.,E",C* 

3065  IF  0*1 1,13  THEN  GOTO  2035 

3066  INPUT  "ENTER  BEGINNING  AND  ENDING  DATA  POINT  NUMBER  3 1  FROM  NE  .  JET  IF  E  RAN 
BED  3  '' ,  St  art ,  F  i  n  i  sh 

3067  RED  I M  UxlF  1  m  sh-St  »r  t  +  1  '• ,  CP  i  m  sh-3*  »r*  ♦!  T 

3068  J-I 

3069  FOR  I »$»  art  TO  Finish 

3070  Ux< J>*X< I > 

3071  Vy  <  J ) »Y< I ) 

3072  MOVE  X  (  I  >  ,  Y  <  I  > 

3073  DRAW  X<U,0 

3074  PENUP 

3075  J«J+1 

3076  NEXT  I 

3077  Eps i I on»IE-6 

3078  CALL  In*  ur.O-1 ,  Epsi  !  on,  U'  •  *  • ,  Vv .  *  • ,  Int  «.;r  si  > 

3079  PRINT  USING  3030;  U*  ■  1  ■ ,  U«  <  J- 1  ■ .  Ir.*  <  •}*•  si 

3080  IMAGE  "THE  AREA  REPRESENTED  BY  THE  SHADED  REGION  BETWEEN  DOMAIN  '  AwUES  " . 
DODD.  SO,  "  AND  "  ,  DDDD.  DD,  "  IS  ,  MD.  DDDE 

3001  PRINT  LlN<2i 

3083  LINPUT  "DO  YOU  WANT  GRAPHICS  FPItlTOUT  VIA  THE  TheF  MAl.  !e I riTE=  ’ “ .  L * 

3090  IF  <L*t  I.  U*“Y”>  AND  <  Cod*  *  I '  THEN  DUMP  GRAPHICS 

3093  EXIT  GRAPHICS 

3100  GCLEAR 

3110  SUBEXIT 

3120  SUBEND 

3130  SUB  Pr  intit  A*,  BS,  Xt»>  ,Y->  ♦>,  N  • 

3140  OPTION  BASE  1 
3130  PRINT  PAGE 

3160  PRINT  "DATA  POINT  RENUMBERING  FOR  BASELINE  CORRECTED  ,  EXPANDED  PLOT" 

3170  PRINT  A* 

3180  PRINT  B* 

3190  PRINT 

3200  PRINT  "DOMAIN  .'ALOE  RANGE  VALUE  NEW  DATA  POINT  NUMBER" 

3210  PRINT 

3220  FOR  I«1  TO  H 

3230  PRINT  USING  2230:  X’- !  ' ,  Y  ■  I  >  .  I 

3240  NEXT  I 

3230  IMAGE  DDDD.  ddxxxxnnnxxxxx.  SD.  BDXXNX:  XXXXXX,  DDD 
3260  PRINT  L I N i 6  ’ 

3270  SUBEXIT 
3230  SUBEND 

3400  SUB  Int  un*,N,  Eps ,  X<  »)  ,  YC  »)  ,  In*  >  1  Cubic  Spline  Integrator 


3430  B»Udt »*<N< «3 •  OP  < X*  1 )  '■X<N!  j  OR  Eps  *0) 

3490  IF  B sdd* *»0  THEN  3600 

3320  PRINT  LIN<.2  ■,  "ERROR  IN  SUBP=OCRPM  Intun,  " 

3330  PRINT  USING  3340; N. Eps . X, l > , X- N  ■ 

3340  IMAGE  "N»" , DDD, SX, "Eps«" , MD. 6DE ,  .  I .  M2 . 6DE  ,  3  ,  "  ■**•*',  *-'D .  6 DE , 
3330  PAUSE 


Jk  • 
% 

>  ■ 


i 
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3560  GOTO  3430 

3600  OPTION  BA3E  t 

3t‘10  DIM  B621N-1  >,S>,rO,G> 2:  N-l  ■ 

3610  FOR  1-2  TO  N-l 
3630  m=13v  I *-1  >-X<  I-I  > 

3640  X*X<I>-X6I-1  > 

3650  B<I>».5*X'H 

3660  T*'  <Y< 1*1  .>-<•'  1)  >■  <X<  1*1  I  .<  - 1 V  ■  I  >-■<•  1-1  •>  >  ■  H 

3670  S<I>-2*T 

3630  G>I)*3*T 

3690  NEXT  I 

3700  9  1 —  3 1  N  1 3 0 

3720  W*S-4*SQP< 3> 

3740  IJ  =  0 

3760  FOR  1*2  TO  N-l 

3770  T*W*< -S< I >-BiI >*$< 1-1 . 5-E< I  •  >*S'  1*1  »G‘  I >  • 

3780  N*ABS‘ T ) 

3600  IF  H,‘U  THEN  U  =  H 

3820  S<I'-*S(I  <*T 

3330  NEXT  I 

3850  IF  U--Epi  THEN  3740 

3860  I nl *3 

3920  FOR  1*1  TO  N-l 

3930  H*X  f I *1 >-X< I  4 

3940  tnt«lnt*<.S»H*<Y-  I  >*Yi  1*1  •  -  1 .  24*H'3»<8-  !'•»$•  I*  l  ■ '•  > 

3950  NEXT  I 
3960  SUEEND 

4030  SUB  Ki  ne  1 1  cs  ’  AS ,  BB,  X<*  >  ,  Y\  r  ,N.Arei-  ’  N-:-r  -  j  c-t  h*rn.  s'  D3C  Ftneticj 

4010  OPTION  BASE  1 

4015  ON  ERROR  GO SUB  Trip 

4020  DIM  LogS'N>,F.-  ict<41N; .  T*i.:p<N),A- 4:N.  ,D'  N  ,DYri:»  5 !  N  ’ .  C  N  ' 

4030  CALL  Pr  in*  I'.fl*,  B#,  Xi  Y'.*  . ,  N) 

4040  PR  I  NT  PAGE 
4050  PRINT  L  I  N  v  2  -* 

4060  PRINT  91 
4070  PRINT  B* 

4030  PFINT  "Kinetics  d*t*rmi  na' -  on" 

4090  PRINT 

4100  PRINT  "  T  €  p '  V  i  1  -■  T  Ar-1  l  Fr  i:  t  .  LN-  ord:  r.  It,  t  '  P«  I  11-  1  ■  '4  P 

it«  P’  i " 

4110  PRINT 

4130  MOT  T«»p*<r>.'X  !  Cel  cut  St «  1  TCIO 

4132  FOR  1*1  TO  N  'INSURE  NON  NEGATIVE  ARGUMENT 

4134  IF  Y  \  I  t  *  1  THEN  GOTO  4140 

4136  Legd'  I '*LOG< V'l  ■  > 

4138  GOTO  4145 
4143  Logd> I >*@ 

4145  NEXT  I 

4180  FOR  1*1  TO  -3  'SET  UP  ARRAYS 
4190  C  Y I  •  »X 1 1 > 

4200  0< I >-Y6 I 

4210  NEXT  I 
4220  Ep*tlon»lE-6 

4230  FOP  1*4  TO  N  'CALCULATE  PARTIAL  AREAS 
4235  RED  IM  C  ••  I  ) ,  D <  I  ' 

4240  C'I)*X<.I> 

4250  DC I  >«Y< I > 

4270  CALL  Int  un<  I ,  Epsi  t  on,  C<  •>,  Dc*  1 ,  P»r  t  i  sl_sr«  i: 

4280  Fnct  <  I  )*Fir*  i  ii_ir«*  flrti 
4285  DISP  Psrtisl  *-•  i.  Fr  sc  1 1 1 '  ,  l 
4*.  90  NEST  I 
4300  FOR  1*5  TO  N 

4  310  3fr«<'I'*'Fr»ctfl  -Fr  ic  * <  1-1  >  •  ’ . :C6  1  , -C  '  1  - 1  ?  - 
4320  NEXT  I 
4330  FOR  1*5  TO  N 

4340  PRINT  USINC  43601  >  ■  I  ■ .  Tf.p  I  >.F- ic»  <.  I  ■  ,Ligd'  I  .  .  Dt>- let  '  I  >,  1 
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4350  NEXT  I 
4355  PRINT 

4360  IMAGE  BDDBXXXX,  D .  DBEXXX,  2.  DB9XXXXXXXXXX>  HD .  DEDXXXXXXXXXXX ,  E .  BZ'EXX  BOB 
4370  !  PLOT  EYract  v»  Fract  AnB  CHOOSE  LINEAR  SECTION  FOR  KINETICS 

4380  PRINT  PAGE 
4330  PRINT  L I N  **  3  > 

4400  PRINT  "KirNtKi  plot  of  dxrdT  vi  x.  x  ■  fraction  d«ccMpoi«d" 

44I0  PRINT  LIN<3) 

4440  C«FMM*xy<Dfr»c» <*>,N> 

4450  D  =  RNMi  nv<Dfract  Of  ,N> 

4400  E*.l 
44/0  F=, C-B 1 ' 10 

4475  INPUT  “ENTER  PLOTTER  CODE:  1=CRT;  2-9S72A" ,Cod» 

4480  IF  Cod«=l  THEN  PLOTTER  IS  “GRAPHICS" 

4435  IF  Cod*»2  THEN  PLOTTER  IS  ?,5."9372A" 

4490  GRAPHICS 

450O  LOCATE  15,90,15,75 

4510  SCALE  0, l.D.C 

4520  AXES  .2,2»F,0,D 

4530  FOR  1=5  TO  N 

4540  MOVE  Fractal, Dfractdl 

4550  LABEL 

4560  NEXT  I 

4578  LORO  6 

4580  FOR  I«@  TO  1  STEP  .2 
4590  MOVE  I,D-.1*F 
4600  LABEL  USING  "Z.BBB“;I 
4610  NEXT  I 
4628  L0RC  8 

4630  FOR  I »D»2»F  TO  C  STEP  2»F 

4640  MOVE  -.01, I 

4650  LABEL  USING  "P.DDE“;I 

4660  NEXT  I 

4670  LORC  6 

4680  MOVE  . 5 , D- . 3*F 

4690  LABEL  "Fraction  D*co»po*«d" 

4700  LttR  PI  -2 

4710  MOVE  .  025,  B  =  5*F 

4720  LABEL  “R»l  at  i  u*  Rat*’ 

4730  LOIR  8 
4740  PAUSE 

4750  IF  Co  d«*  I  THEN  BUMP  GRAPHICS 

4760  CCLEAR 

4770  EXIT  GRAPHICS 

4780  PPINT  PAGE 

4790  PRINT  LIN<2> 

4800  PRINT  "LINEAR  REGRESSION  ANALYSIS  OF  BATA" 

4818  PRINT  L I N r 2 1 

4820  INPUT  “ENTER  BEGINNING  ANB  ENDING  POINTS  FOR  CALCULATION" . B«gi n. End 

4030  Rt*»£nd-B«girvM 

4040  REDIM  T2 1 Rt * ) , R2 1 Rt a > 

4850  IMAGE  “Plot  i»  b«tw*«n  tt»p«ratur*a  ",BBDB."  and  “,DDBB,“  Ktlvin" 

4868  FOR  I* l  TO  Pt» 

4861  T2'  I  i«T»»c  Bfj1n*l-1) 

4862  F2d  >»Logd<  B«gi  n*I-l  f 
4965  NEXT  I 

4866  PRINT  USING  4850* l'T2v l >, lrT2<Pl»> 

4870  PRINT 

4880  CALL  Lir>«ar<T2C*).F2r*>,Rt»,A,B,C0,Cl,C2.C3.C4,C5,C6,C7,ce.C9,V> 

4890  PRINT  USING  4900; -B»l .  9872  >  Act  i  vat  i  on  inr  gv  in  :»I.r»ol« 

4395  PRINT  USING  4905JV 

4988  IMAGE  "Calculand  acti-'aticn  *n«rgv  tllnaar  r»gr«»»ion)  ia  ",B,BBE,  '  ca'  N 
o  I «  " 

4905  IMAGE  "M«an  »lop«  oarianc*  ii  ”,2.DBDE,"  cal  *ol*“ 

4910  PRINT  LINOj 

4920  ::nia**FNMax2'  tj.  *  >,  i  ,p,  *> 
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493$  :»i  n*FNMi  nZ(  T2'  *'* «  1  ,  P*  4  ' 

4940  /«.*::*FNM*-2<F2<»'-,  1  ,P‘*> 

4950  Yn.in»FNHtr>2<F2<*>,  t .  R*  *  ' 

4960  £*'  X«.»-,-X*tn>'  10 

4970  F*'.  Y»»*-Y«m  n>  '  10 

49S0  IF  Cod«*l  THEN  PLOTTER  IS  "GRAPHICS" 

4985  IF  Cod»*2  THEN  PLOTTER  IS  7,5, "98728“ 

4990  GRAPHICS 

5000  _0C ATE  15.90,15,75 

5010  SCALE  X«.ir.IX»»x>  Vmin.VBi*. 

5020  AXES  24E.  2*F,  X»i  n,  Ytui  n 
5030  FOP  1*1  TO  Pt» 

5040  MOVE  T2< I;,F2<1> 

5050  LABEL  *♦“ 

5080  NEXT  I 

5070  FOP  J*1  TO  Pt* 

5080  F2*  J ' *5*  T2<  J>  +  A 

5090  NEXT  J 

5895  MOVE  T2<  1  ‘,F2<.1> 

5100  FOR  J*2  TO  PtJ 
5110  DRAW  T2<J  >,F2<3/ 

5120  NEXT  J 
5140  LORG  6 

5158  FOR  I*Xmin  TO  STEP  2*E 

5160  MOVE  I,  Y*ir>-.  t»F 
5170  LABEL  USING  "D. DDD" ; 1 «1O00 
5180  NEXT  I 
5200  LORG  3 

5210  FOR  I*V«in+2*F  TO  Y»*.<  STEP  2*F 
5220  MOVE  X*tn-.04*E,  t 
5230  LABEL  USING  "MD.DDD";I 
5240  NEXT  I 


LORG  6 

MOVE  Xi»irt*5*E,Vmo-F 
LABEL  "1  T.tO  v  1000" 

MOVE  Xm  n*.  4*E.  V»i  n*5*F 
-DIF  PI -2 

„  ABEL  "  1  n  <di  »pl  K«»»nt  >" 

R  AUSE 

IF  Cod«*l  THEN  dump  graphics 
gclear 

EXIT  GRAPHICS  „„  „ 

L INPUT  "DO  YOU  WANT  TO  REPEAT  THE  CALCULATION0", Repeat* 

IF  R«p«»t*Cl,  n*"V*  THEN  GOTO  4760 
SUBEXIT 

Tr  *p:  DISP  ERRM* 

PAUSE 
GOTO  4040 

p  DEF  FNM»xv<XT*  ',N' 

I  Hoxv*X<5> 

I  FOR  1*5  TO  N 

I  IP  X'S)>H«<y  THEN  M*-!V*X'.I  > 

)  NEXT  I 
»  RETURN  M»xy 
)  FNENS 

1  DEF  FNMioy<X<*),N> 

1  Mi  nv"Xi'5 ' 
i  FOR  1*5  TO  H 

1  IF  XCIXMiny  THEN  Minv*X'l' 

»  NEXT  I 
>  RETURN  M«ny 
9  FNCND 

l  SU»*u^«,r'X<»'-.Y'*:  ,N, A. B,R«94».R«»t». Tot »l*l,R«0P»*,R«»»*,F,DFr«9. I»f 
,  Abort  ,  V*r  i  *rc  t  • 
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9710  '  ♦*»  MODEL:  Y»A*B»X  *** 

57*0  ON  ERROR  GOTO  Bctnb 
5700  Abor*  «0 
5740  Y1*X1*Z*X2“Y2»6 
5750  FOP  1*1  TO  N 

5700  X1*X1+X<  I  >  [  Sum  of  X  * 

5770  V  l  *  Y 1  ♦  Y  t  1  >  I  Sum  of  Y  s 

57S0  X2*X2*X< I >*Xt I >  <  Sum  of  X  tqutrtt 

57S0  Y**Y2+Y< l >*Y< I >  !  Sum  of  Y  square* 

5S00  :«2*Xtl)*m)  !  Sum  of  XY's 

5810  NEXT  I 

5820  Xl-Xl'N  !  Mean  X 

5800  Y1»Y1  N  •  Mean  Y 

5840  B»(Z-N«X1*Y1  >-'<X2-N*Xl*Xl  ' 

5350  A*Y1-B*XI 

5880  Total  ss«Y2-N»Yl*Yt  •  Total  Sum  of  Squares 

5370  Regu  ••;-N*Xl*Yl  >"2  '<X2-N»:;i*Xl ) !  Regression  Sum  of  Squares 
5830  PessseTotalss-Regss  <  Residual  Sum  of  Squares 

5830  Regms*Pegss  !  Regression  Mean  Squares 

500O  ResMS»Ressi  <N-2>  !  Residual  Mean  Squares 

5005  Var i anc e »SQR<Res»S''<X2-N*XleXl > >  I  Slope  Variance 
5910  OFF  ERROR 
5920  DEFAULT  ON 

5930  F«R*gMS'Rese.s  •  F  Ratio 

5*40  DEFAULT  OFF 

5950  Dfreg«l  1  Degrees  of  Freedoa 

5980  Dfres*N-2 

9970  Dftot *Dfreg*Dfres 

5930  SUBEXIT 

5990  Bomb:  Abort «1 

8000  SUBEND 

€010  DEF  FNM 1*2 'Xtei, A, B> 

8020  Ma<2*X(A> 

8030  FOR  I*A»1  TO  B 

6040  IF  xav  MaxZ  THEN  Hax2*X<I> 

6050  NEXT  I 
6080  RETURN  Ma>  2 
60~0  FNEND 

6080  DEF  FNMin2tX<ei,A, B> 

6090  Min2»X’A> 

6100  FOR  !•«♦!  TO  B 

6110  IF  X<IXMin2  THEN  Mm2*X<I> 

6120  NEXT  I 
6130  RETURN  Mm2 
61*0  FNEND 

6200  SUB  Pol  yno»  ial<Xfe),Y(*),N,  Degree  ,Coeff*t*'',Reg*s,Res**.Toialss,  Reg** .  Re  sms 
, F,D‘"  re  g.Dfres, Dftot,  Abort  ) 

6210  ON  ERROR  GOTO  BoMb 

6220  Abort »0 

6230  OPTION  BASE  0 

6248  DIM  M*tr i \< Degree , Degree), InuCDegree, Degree' , BfDegree) 

6250  PEDIM  Coeff st Degree > 

62€0  IF  Degree.N-2  THEN  SUBEXIT  ‘  Check  for  higher  degree  than  possible 
6270  Dfreg*Degree 
6280  Df res«N-t-Degr*e 
6290  Dftot»Dfreg*Df res 

6300  FOP  K»0  TO  Degree  1  Set  up  systea  of  equations 

6310  FOR  J«N  TO  Degree 
6320  Matr I x<K,  J)*0 

6830  FOR  1*1  TO  N 

8  340  Matrix’ K, J>-Matr i xtK, J»*FNGvK>eFHG<  J> 

6250  NEXT  I 

6360  Ma’ri  <t J, K >*Nat r i xf PC , J) 

6370  NEXT  J 

63  80  BflO*0 

8390  FOR  1*1  TO  N 
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6400  BOO-BOO+YCI  '*FNC'L  • 

6410  NEXT  I 
6420  NEXT  K 

6430  MAT  Ir.u*INV>;i1*tr  j  ;;> 

6440  MAT  Coeffi*Inw*B 
6450  FOP  I«1  TO  N 
6460  X1*X1«X<I> 

6470  X2»X2+X<I>*X<I> 

6430  Yt»YI*V<I> 

6400  Y2»Y2*Y<I)4Y(I) 

65O0  Z«2*X< I )*Y< 1 ) 

6510  NEXT  I 
6520  Yt«Yl.'N 
6530  Xl-Xl- N 
6540  Tot»U*«Y2-N*Yl«Yl 
6550  GOSUB  Regst 
6560  Fes*j«Tota! ss-Regss 
6570  Peg*s«R«gss-  Bfreg 
6530  Res*s*Re sss ' Df rej 
6590  OFF  ERROR 
6600  DEFAULT  ON 
6618  F»Reg»s  Resas 
6620  DEFAULT  OFF 
6630  SUBEXIT 
6640  Regss:  Regss*0 
6650  FOR  I»1  TO  N 
6660  J-0 

6670  FOR  L«0  TO  Degree 

6680  J«J*X<I)"L»CoeffsCL> 

6690  NEXT  L 

6780  Regs*«Regs*+<.  J-Yl  >*2 

6710  NEXT  I 

6720  RETURN 

6730  SUBEXIT 

6740  DEF  FNG<m;>«X<.I>“N 

6750  Bomb:  Abort »1 

6760  SUBEND 

6770  SUB  Quad<X<e),Y«:«),Z) 

6780  PRINT  LIN<2> 

6785  OPTION  BASE  1 
6790  DIN  A<3> 

6800  CALL  Polynonl »1 <X<»> , Y< , Z, 2,  AC *> , CO, C 1 , C2, C3, C4, C5, C6, C7, C8, C9' 

6310  PRINT  USING  6820; -A<2>*1 . 9872, A< 1 >* l . 9872 

6311  PRINT  A<3> 

6020  IMAGE  "THE  QUADRATIC  COEFFICIENTS  ARE  "  ,  MD.  DDE,  *  and  *  ,  ND.  DDE 
6830  PRINT 
6040  SUBEXIT 
6050  SUBEND 


1  Soloe  the  i;atii)  of  equations 


1  Tot *1  Sun*  of  Square! 

!  Regression  Sun  of  Squtrti 
!  Residual  Sun  of  Squares 
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NOMENCLATURE 

A  Total  Area 

a  Partial  Area 

D(t)  Ordinate  Deflection 

E*  Activation  Energy 

f(x)  Algebraic  Function  of  Fraction  Reaction 
G  Gibbs  Free  Energy 

H  Enthalpy 

K(T)  Equilibrium  Constant 

k(T)  Rate  Constant 

R  Universal  Gas  Constant 

T  Temperature 

Tm  Temperature  of  Maximum  Exotherm 

t  Time 

X  Fraction  Reacted 

Z  Preexponential  Factor  in  Arrhenius  Equation 
Heating  Rate 
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